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Foreword 

This  report  describes  results  obtained  during  a  three  year  research  study  on  nonlinear  mod¬ 
eling  and  control  of  flexible  space  structures  with  application  to  rapid  slewing  and  precision 
pointing  of  space-based,  directed  energy  weapons.  The  project  is  funded  by  SDIO/IST  and 
managed  by  AFOSR/SDIO  (AFSC).  Results  reported  herein  are  for  the  period  1  Sept.  1987 
-  31  Aug.  1990.  The  research  effort  was  conducted  at  Techno-Sciences,  Inc.  Greenbelt, 
Maryland  office.  The  project  principal  investigator  was  Dr.  William  H.  Bennett,  Director 
of  Control  Systems  Projects  at  TSI.  Coinvestigators  included  Drs.  Gilmer  Blankenship  and 
Harry  Kwatny.  Dr.  Kwatny  is  Vice  President  and  member  of  Board  of  Directors  at  TSI 
and  is  Raynes  Professor  of  Mechanical  Engineering  and  Mechanics  at  Drexel  University.  Dr. 
Blankenship  is  Vice  President  and  member  of  Board  of  Directors  at  TSI  and  Professor  of 
Electrical  Engineering  at  University  of  Maryland,  College  Park,  MD.  Associate  investigator 
on  the  project  was  Dr.  Oussima  Akhrif,  who  is  currently  parttime  member  of  TSI  control 
systems  staff  and  visiting  assistant  professor  at  Case  Western  Reserve  University,  Cleveland, 
Ohio. 

During  the  first  two  years  the  project  was  managed  by  Ltn.  Colonel  J.  Crowley/ AFOSR 
and  Dr.  A.  Amos/AFOSR.  We  wish  to  thank  both  of  these  individuals  for  their  insight  and 
direction  on  this  project.  During  the  last  year  the  project  technical  manager  was  Dr.  Alok 
Das/AFAL. 
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1  Research  Objectives  and  Project  Summary 

This  report  describes  results  of  a  three  year  project  on  Nonlinear  Dynamics  and  Control  of 
Flexible  Structures  whose  principal  objective  was  to  develop  modeling  and  control  design 
methods  for  nonlinear,  distributed  parameter  dynamics  arising  in  requirements  for  rapid 
slewing  and  precision  pointing  of  multibody  space  structures  with  elastic  interactions.  The 
modeling  effort  centered  around  the  development  of  a  generic  class  of  models  for  the  struc¬ 
tural  interactions  effecting  laser  system  line-of-sight  (LOS)  pointing  for  a  Space-Based  Laser 
(SBL)  weapon.  Control  studies  focused  on  the  problem  of  rapid,  multiaxis  large  angle  slew¬ 
ing  of  the  generic  SBL  system  with  requirements  for  precision  system  LOS  pointing.  One 
objective  for  this  project  was  to  establish  a  framework  for  experimental  validation  of  the 
class  of  nonlinear  control  laws  under  consideration.  To  the  extent  possible,  simulation  mod¬ 
els  were  based  on  predicted  dynamics  and  control  configurations  for  the  first  test  article 
currenlty  under  construction  at  the  ASTREX  facility  at  Air  Force  Astronautics  Laboratory. 
Simulation  results  are  included  for  several  possible  experiments.  A  proposed  experimental 
protocol  for  validation  of  compensation  of  nonlinear  couplings  effecting  slewing  and  precision 
pointing  is  briefly  described. 

Our  approach  for  slewing  and  pointing  control  was  to  develop  several  practical  extensions 
and  applications  of  the  theory  of  nonlinear  control  system  synthesis  based  on  linearization  by 
feedback  transformation.  These  methods  include  control  laws  based  on  “dynamic  inversion” 
such  as  the  so-called  “method  of  computed  torques”  popularized  in  robotics.  However,  our 
analysis  indicates  that  certain  essential  limitations  arise  in  control  of  distributed  parameter 
dynamics  arising  in  multibody  systems  with  elastic  structural  interactions.  The  approach 
we  consider  is  based  on  the  implementation  of  a  dynamic  inverse  of  the  response  of  certain 
principal  system  outpi  "s  representing  multibody  system  pointing  and  alignment  to  system 
actuation.  As  will  be  discussed  in  the  sequel  this  concept  of  input/output  linearization  (when 
it  can  be  achieved)  involves  decoupling  certain  complementary  or  “zero  dynamics”  from  the 
principal  system  outputs.  Since  the  decoupled  dynamics  will  in  general  be  nonlinear  we  refer 
to  this  class  of  control  laws  as  Partial  Linearizing  Feedback  (PLF). 

Several  critical  system  requirements  and  control  systems  can  be  identified  for  a  prototype 
SBL  system.  Precision  line-of-sight  (LOS)  pointing  and  tracking  must  be  achieved  in  concert 
with  requirements  for  rapid  slewing  of  spacecraft  primary  body  and/or  articulation  of  weapon 
system  aperture  (see  Figure  1.1).  Optical  beam  quality  and  jitter  performance  will  be  limited 
by  extent  to  which  precision  multibody  alignment  of  the  structure  supporting  the  optical 
train  can  be  dynamically  controlled  during  and  after  spacecraft  system  slewing  maneuvers. 
Thus  slewing  and  pointing  control  for  an  SBL  system  will  involve  the  integration  of  several 
critical  control  systems.  An  important  feature  obtained  by  PLF  control  is  decoupling  of 
control  loops  for  processing  slewing  and  pointing  commands.  The  class  of  PLF  control  laws 
considered  in  this  study  for  LOS  slewing  and  pointing  admits  a  special  structure  which 
permits  the  required  decoupling/linearizing  transformations  to  readiliy  implemented  with 
multiple  actuators.  One  way  in  which  this  may  be  used  to  advantage  is  in  the  integration  of 
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SPACE  BASED  LASER  POINTING  AND  TRACKING 


Figure  1.1:  Generic  Space-Based  Laser  System  with  Multiple  Control  Systems 


continuous  mode  (e.g.  reaction  wheels,  CMG’s,  etc.)  and  discontinuous  (on/off)  actuation 
(e  g-  jets). 

Outline  of  the  Project.  The  first  year  effort  included  a  survey  of  the  available  literature 
on  the  dynamic  models  and  control  problems  for  SBL  systems.  A  generic  class  of  SBL 
models  was  developed  including  multibody  dynamics  of  a  dual  mirror  laser  beam  expander 
with  structural  flexure  arising  in  the  metering  truss  which  supports  the  secondary  mirror 
relative  to  the  primary  mirror.  Initial  control  studies  focused  on  large  angle,  but  single  axis 
slewing  maneuvers  and  restricted  the  model  to  planar  motions.  Computer  simulations  were 
completed  demonstrating  the  slew  rates  and  control  requirements  for  PLF  control.  During 
the  second  year  we  expanded  the  computer  simulation  model  to  include  multiaxis  slewing 
motions  and  demonstrated  several  alternative  control  laws  based  on  PLF  compensation  which 
included  both  continuous  and  discrete  actuation  modes.  A  goal  of  the  effort  was  to  establish 
a  basis  for  experimental  validation  of  the  slewing  control  laws.  Extensive  tradeoff  studies 
were  performed  with  simulation  models  scaled  to  represent  predicted  dynamic  response  for 
the  first  test  article  currently  under  construction  at  the  ASTREX  test  facility  at  Air  Force 
Astronautics  Laboratory. 

New  Results  Obtained.  A  potentially  important  feature  of  feedback  linearization  is  the 
extent  to  which  the  idea  can  be  integrated  with  other  standard  design  methods  for  multivari¬ 
able  control  systems.  We  have  described  several  alternatives  which  indicate  advantages  of 
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the  approach  for  control  of  flexible  structures.  In  particular,  the  integration  of  active/passive 
techniques  for  vibration  suppression  can  be  readily  incorporated.  The  idea  of  deformation 
shaping  control  has  been  developed  by  Dr.  T.A.W.  Dwyer  as  part  of  this  project.  The 
method  is  reported  in  several  papers  prepared  as  part  of  this  project.  PLF  compensation  of¬ 
fers  a  framework  for  resolving  certain  design  tradeoffs  for  the  use  of  active/ passive  vibration 
control  for  rapid  slewing  and  precision  pointing  of  SBL  (or  other  multibody  systems). 

Robust  stabilization  and  control  performance  of  nonlinear  systems  is  an  important  issue 
in  many  applications.  For  the  current  study  we  are  interested  in  design  of  control  laws 
whose  stability  and  performance  can  be  predicted  and  validated  in  experiments.  Several 
techniques  for  robust  stabilization  of  nonlinear  systems  have  been  developed  by  Spong  and 
Vidyassagor,  Slotine  and  Sastry,  Corless,  Gutman,  and  others.  These  methods  all  rely  on 
special  structure  of  the  dynamic  model  uncertainty  called  “matching  conditions”.  In  this 
research  we  have  developed  a  new  approach  which  utilizes  the  framework  of  adaptive  control 
to  provide  robust  stabilization  of  nonlinear  systems.  The  approach  does  not  rely  on  structure 
matching  conditions  and  represents  a  significant  improvement  over  available  methods. 

Overview  of  the  Report  for  Third  Year  Activities.  In  this  report  we  start  with  a 
review  of  the  basis  for  nonlinear  PLF  control  design  with  emphasis  on  stabilization.  We 
demonstrate  several  alternative  implementation  schemes  for  PLF  control  which  can  arise 
from  several  different  problems.  In  particular,  we  show  how  PLF  compensation  can  be  ob¬ 
tained  using  switching  mode  actuation  by  design  of  sliding  mode  control  systems.  We  also 
indicate  how  techniques  from  adaptive  control  can  be  used  to  improve  robustness  when  dy¬ 
namic  models  are  subject  to  parametric  uncertainty.  We  next  develop  a  complete  framework 
for  modeling  multibody  systems  with  elastic  structural  interactions  using  Lagrange’s  equa¬ 
tions  for  continuum  mechanical  systems.  We  specialize  these  systems  to  the  generic  SBL 
models.  We  then  discuss  a  framework  for  computation  of  PLF  control  laws  for  Lagrangian 
systems  which  focuses  attention  on  the  notion  of  nonlinear  system  zero  dynamics.  For  prac¬ 
tical  design  of  control  for  large  space  structures  we  describe  implementation  of  PLF  controls 
for  reduced  order  models.  We  specialize  the  standard  model  reduction  problem  for  PLF 
control  by  considering  time  scale  separation  of  the  system  zero  dynamics.  Finally,  we  briefly 
describe  a  protocol  for  experiments  which  can  be  used  to  validate  PLF  control  performance 
predictions.  Simulation  results  are  included. 

The  first  [BBKA88]  and  second  [BBKA89]  year  reports  contain  additional  results  on 
control  design  and  tradeoff  studies.  '■ 

Professional  Personnel.  The  principal  investigator  for  this  project  is  Dr.  William  H. 
Bennett  and  co-principal  investigators  are  Drs.  Harry  G.  Kwatny  and  Gilmer  L.  Blankenship 
from  TSI.  Dr.  Thomas  A.  W.  Dwyer  was  consultant  on  the  project.  We  would  also  like  to 
acknowledge  the  parttime  support  from  Dr.  Oussima  Akhrif  who  graduated  in  July  1989 
with  Ph.  D.  from  Electrical  Engineering  Department  of  the  University  of  Maryland.  Dr. 
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Akhrif’s  dissertation  developed  aspects  of  feedback  linearization  for  the  multibody  SBL 
models  developed  in  this  study  and  some  important  results  are  summarized  in  this  report. 

Technical  Reports/Presentations.  As  part  of  the  research  program  we  have  organized 
two  invited  technical  sessions,  presented  several  technical  papers,  and  submitted  several 
additional  papers  for  publication  as  follows. 

1.  Nonlinear  Dynamics  and  Control  of  Aerospace  Systems,  invited  session  at  27th  IEEE 
Cntrl.  Dec.  Conf.,  Austin,  TX,  Dec.  1989. 

2.  Robust  Control  of  Uncertain  Nonlinear  Systems,  invited  session  at  1989  Amer.  Cntrl. 
Conf.,  Pittsburgh,  PA,  June  1989. 

Publications/Presentations 
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2  Feedback  Linearization  and  Stabilization  of  Nonlinear  Systems 

Conventional  techniques  for  stabilization  of  nonlinear  systems  via  feedback  control  are  still 
very  limited  and  tend  to  be  tailored  to  specific  situations.  Among  the  most  promising  new, 
general  approaches  utilize  linearization  (local  or  possibly  globed)  by  Exact  Feedback  Lin¬ 
earization  ( EFL )  [HSM83,  KC87],  EFL  methods  are  based  on  earlier  work  of  Krener  [Kre73] 
and  Brockett  [Bro7bj  which  demonstrated  that  a  large  class  of  nonlinear  dynamical  systems 
can  be  exactly  (i.e.  globally  linearized)  by  a  combination  of  nonlinear  transformation  of 
the  state  coordinates  with  nonlinear  state  feedback.  More  recently,  the  connection  between 
these  methods  and  the  idea  of  input-output  (or  Partial)  Linearizing  Feedback  (PLF)  by  con 
struction  of  a  system  inverse  [Hir79]  has  been  articulated  in  a  series  of  papers  by  Byrnes 
and  Isidori  [BI85,  BI84].  These  connections  have  engendered  a  series  of  design  methods  with 
representative  -esults  for  specific  applications  by  Kravaris  and  Chung  [KC87]  and  Fernan¬ 
dez  and  Hedrick  [FH87].  In  this  section  we  will  show  how  fundamental  these  construction 
can  become  in  control  system  design,  discuss  alternatives  for  implementation,  and  suggest 
some  approaches  to  integrating  the  nonlinear  design  philosophy  with  more  conventional  ap¬ 
proaches.  We  focus  attention  in  this  section  on  fundamental  concepts  culminating  in  the 
description  of  the  design  approach  for  multibody  systems  from  the  perspective  of  Lagranian 
mechanics. 

The  idea  behind  feedback  linearization  is  conceptually  simple.  We  start  with  a  nonlinear 
system  model, 

x  =  f(x)  +  G(x)u,  (2.1) 

y  =  M*)»  (2-2) 

where  x  6  3?n,  u,y  6  with  G  =  [<7i,. . .  ,gm ]  and  assuming  the  vector  fields  /,  <7,  are  C°° 
for  each  i  —  1, . . .  ,m  and  /( 0)  =  0.  The  model  structure  assumes  that  the  control  u  enters 
linearly.  The  feedback  linearization  problem  is  to  find  a  change  of  basis  in  the  state  space, 
z  =  T(x),  with  T  diffeomorphic  and  a  feedback  law, 

u  =  a(i)  -f  0(x)v, 

such  that  in  the  new  ( z,v )  coordinates  the  (closed  loop)  model  has  the  form, 

z  =  Az  -f  Bv. 

We  remark  that  if  it  possible  to  find  such  a  control  law  then  the  linearization  is  achieved 
thr  ough  the  introduction  of  active  control  authority.  An  important  feature  for  control  system 
design  is  that  the  range  of  validity  of  the  linearization  is  given  by  the  transformations, 
T(x),a(x), /3(x).  The  functions  may  be  defined  locally  or  globally. 

In  contrast,  the  conventional  approach  to  control  design  would  be  based  on  a  linear  model 
obtained  by  Taylor  expansion  of  the  vector  fields  about  given  equilibrium  conditions ;  xtg,ueq, 
satisfying, 

0  =  f(xeq)+  G(xeg)ueg. 
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The  conventional  linear  model  used  for  control  design  represents  perturbation  dynamics  with 
respect  to  the  equilibrium  conditions  and  assumes  the  form, 

(Ax)  =  A  Ax  -f  B  Au, 

y  =  C  Ax 

where  Ax  =  x  —  xe<r,  Au  =  u  —  ueq  and 


.  B  =  G(z)U„  .  C 

X  —  Xfq 


8h 

dx 


z=*e. 


A  major  source  of  model  uncertainty  in  linear  control  system  design  arises  from  assumptions 
leading  to  the  linear  perturbation  model.  In  many  cases  it  may  be  difficult  to  estimate  the 
domain  of  attraction  for  the  equilibria.  Indeed,  in  aerospace  applications  the  control  design  is 
often  based  on  a  combination  of  gain  scheduling  to  take  into  account  the  dependence  of  linear 
perturbation  models  on  operating  point  conditions  which  are  subject  to  variation  (e.g.  trim 
conditions  in  aircraft  flight  control.)  For  example,  the  function  of  a  conventional  autopilot 
for  aircraft  is  to  compensate  for  changes  in  trim  conditions  and  provide  stabilization  so  that 
the  pilot  “feels”  a  standard,  linear  response  to  stick  commands. 

Although  the  concept  of  feedback  linearization  in  control  system  design  is  potentially 
revolutionary,  its  application  has  many  antecedents  in  applications.  The  significance  for 
nonlinear  control  of  flexible  space  structures  is  the  emerging  technology  for  active  control 
and  sensing,  the  dynamics  associated  with  the  CSI  technology,  and  the  ability  of  a  com¬ 
prehensive  approach  to  nonlinear  dynamic  modeling  and  control  design  offered  by  the  ap¬ 
proaches  discussed  in  this  report.  Feedback  linearization  functions  in  certain  applications 
in  a  manner  similar  to  gain  scheduling  [MC80],  however,  linearization  is  achieved  about  a 
“nominal  model”  rather  than  about  an  operating  point.  Thus  equilibria  conditions  do  not 
arise  explicitly  in  linearization.  One  view  of  such  a  controller  structure  is  illustrated  in  block 
diagram  of  Figure  2.1.  The  process  linearization  which  facilitates  the  design  of  the  Unear 
controller  is  obtained  by  the  introduction  of  an  Inverse  Force  Model  (IFM)  for  the  nonfinear 
multibody  system.  The  inverse  force  model  transforms  commanded  accelerations,  ae,  into 
equivalent  system  generafized  forces,  /.  Thus  the  Unear  controller  is  designed  to  yield  desired 
system  accelerations  given  the  generalized  coordinates,  q,  and  their  rates,  q. 

Precursors  to  the  idea  of  feedback  Unearization  is  pervasive  in  control  applications  abound. 
With  the  development  of  the  geometric  theory  of  nonfinear  systems,  computational  tools  and 
design  methods  are  becoming  available  to  address  control  system  design  on  a  much  larger 
scope.  It  is  clear  that  the  concept  of  feedback  Unearization  is  pervasive  in  many  fundamental 
control  methods.  Our  study  has  focused  on  the  considerations  for  practical  implementation 
of  feedback  Unearization  for  rapid  slewing  and  precision  pointing  of  aerospace  systems  with 
multibody  and  elastic  interactions.  Implementation  of  feedback  Unearization  definitely  re¬ 
quires  enhanced  control  authority  and  issues  related  to  technology  for  control  actuation  will 
either  enable  or  restrict  its  appUcation. 
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Figure  2.1:  Nonlinear  Control  Concept  Using  Dynamical  System  Inverse 

It  is  often  suggested  that  controllers  designed  based  on  feedback  linearization  may  be 
sensitive  to  model  assumptions  since  the  linearization  is  achieved  by  cancelation  of  cer¬ 
tain  nonlinear  terms  in  the  system  model.  Our  study  of  rapid  slewing  control  of  a  generic 
space- based  laser  model  has  shown  that  modeling  sensitivity  and  robustness  can  be  obtained 
through  judicious  application  of  control  authority.  A  central  issue  in  nonlinear  control  design 
is  the  limits  of  available  control  authority.  For  example,  actuator  saturation  contributes  to 
limits  on  control  authority.  Discontinuous  or  saturation  mode  operation  of  control  actua¬ 
tors  can  often  be  desirable  but  such  considerations  are  often  not  addressed  by  linear  design 
methods.  We  have  shown  that  specific  consideration  for  coordination  of  discontinuous  and 
continuous  modes  of  control  actuation  in  slewing  control  methods  for  multibody  systems 
can  be  readily  found.  Finally,  it  is  clear  that  the  cost  of  feedback  linearization,  in  terms 
of  increased  control  authority,  sensor  measurement  complexity,  or  computational  burden  for 
online  implementation  may  not  always  be  necessary  to  achieve  system  performance  objec¬ 
tives.  We  have  also  demonstrated  that  such  methods  can  be  readily  integrated  with  standard 
approaches  for  linear  design  once  a  primary  system  control  objective  is  identified  in  terms 
of  a  primary  system  output.  As  we  shall  show  the  role  of  conventional  linearization  and 
control  design  can  be  relegated  to  a  subsystem  whose  dynamics  are  decoupled  (by  the  action 
of  PLF)  from  a  set  of  system  primary  outputs. 

2.1  Computation  of  Partial  Linearizing  Feedback  Compensation 

Partial  linearization  derives  directly  from  the  Byrnes-Isidori  normal  form  for  nonlinear  sys¬ 
tems.  The  essentials  of  the  approach  are  most  easily  developed  for  single-input,  single  output 
systems  and  we  will  present  the  approach  in  that  context.  The  theory  for  extending  these 
results  for  multi-input,  multi-output  problems  is  now  complete  and  references  are  included. 

Consider  a  nonlinear  dynamical  system  in  the  form, 

*  =  f{x)  +  g(x)u 
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y  =  h(x )  (2.4) 

where  /,  g  are  smooth  C°°  vector  fields  on  9?n  and  h  is  a  smooth  function  mapping  5?n  — ♦  9i. 
Now  if  we  differentiate  (2.4)  we  obtain 


V  =  +  9ix)u)-  (2-5) 

In  the  case  that  the  scalar  coefficient  of  u  (viz.  § ^g[x))  is  zero  we  can  differentiate  again  until 
a  nonzero  control  coefficient  appears.  The  number  of  required  differentiations  is  fundamental 
system  invariant  which  plays  a  role  in  constructing  a  system  inverse  and  therefore  in  PLF. 
The  Byrnes-Isidori  analysis  shows  that  this  integer  number  is  analogous  to  the  relative  degree 
for  a  linear  system  [BI84]. 

The  above  construction  can  be  made  precise  using  the  notation  of  differential  geometry 
which  has  found  application  in  analytical  mechanics  [Arn78].  We  will  need  only  the  notion 
of  Lie  derivative  and  Lie  bracket.  The  Lie  (directional)  derivative  of  the  scalar  function  h 
with  respect  to  the  vector  field  /  is 

L,(h)  =  l,dh,f)-.=  ^f(x).  (2.6) 

Since  the  above  operation  results  in  a  scalar  function  on  9?n,  higher  order  derivatives  can  be 
successively  defined 

L)(h)  =  L,{L)-\h))  :=  (iL)~'(h),  /).  (2.7) 

Then  we  can  write  (2.5)  as 

y  =  {dh,f)  +  (dh,g)u 

=  Lf(h )  -I-  Lg(h)u.  (2.8) 

If  Lg(h)  =  0  then  we  differentiate  again  to  obtain 

y  =  { dLf(h),f)  +  {dLf(h),g)u 

=  L){h)  +  Lg{Lf{h))u.  (2.9) 

If  2/fl(Z/5_1(h))  =  0  for  k  =  1, . . .  ,r  —  1,  but  Lg(Lrfl(h))  ^  0  then  the  process  terminates 
with 

-jf  =  L’,W  +  L,(L'f'(h))u.  (2.10) 

The  system  (2.10)  can  be  effectively  inverted  by  introducing  a  feedback  transformation  of 
the  form 

<2U) 

which  results  in  an  input-output  response  from  v  — ►  y  given  by 

dry 
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a  linear  system. 

The  integer  r  >  0  can  be  viewed  as  a  relative  degree  for  the  nonlinear  system  (2.3)-(2.4). 
Note  that  if  we  define  new  state  coordinates  z  £  5?r  as 


zk=L)-'(h ),  *  = 


(2.12) 


for  the  r-dimensional  nonlinear  system  (2.10),  then  the  system  model  can  be  written  in  state 
space  form  as, 


z  = 


0 

'  0 

1 

0  • 

•  0 

0 

0 

0 

1  • 

•  0 

0 

0 

0  • 

.  1 

2  + 

• 

0 

0 

0  • 

•  0 

0 

a(*)  +  p(x)u 

(2.13) 


where 


a(x)  =  L’,(h),  p(x)  =  L,(L'f\h)). 


(2.14) 


More  generally,  using  the  new  coordinates  z  (2.12)  and  introducing  a  nonlinear  feedback 
control  of  the  form 

(2.15) 


u  — 


where 


/>(*) 

<-(*)  =  J>i*(t)  +  W, 


r— 1 


(2.16) 


k=0 


p(x)=L,(L)-'(h)),  (2.17) 

with  ft  for  k  =  0, . . .  ,r  —  1  real  positive  coefficients  then  the  equations  (2.8)-(2.9)  can  be 
written  in  ‘reduced’  form; 


'  0 

l 

0  ••• 

0  ■ 

( 0  \ 

0 

0 

1 

0 

0 

z  = 

0 

0 

o 

0 

Z  + 

: 

.  -ft 

-ft 

-ft  ... 

-ft- 1 . 

l  1  / 

y  =  [1,0, ...  ,0]  z. 

2.2  Nonlinear  System  Transmission  Zeros 


(2.18) 

(2.19) 


Note  that  the  process  leading  to  (2.18)-(2.19)  provides  an  equivalent  state  space  realization 
for  the  vi -*  y  input-output  response  of  McMillan  degree  r  <  n  (the  dimension  of  the  original 
state  space  model  (2.3)-(2.4))  by  decoupling  a  portion  of  the  system  dynamics  from  the  output 
response.  This  is  depicted  in  Figure  2.3.  Thus  the  new  state  coordinates  z  are  a  ‘partial’  state 
for  the  system.  Thus  stabilization  of  (2.18)-(2.19)  cannot  guarantee  stabilization  of  the  full 
state  model  (2.3)-(2.4).  We  remark  that  in  the  case  that  h(x)  is  such  that  the  relative  degree 
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r  =  n  then  the  state  space  transformation  (2.12)  for  k  =  1, . . . ,  r  together  with  the  feedback 
transformation  (2.15)  exactly  linearizes  the  full  system  state  model  (2.3).  The  methods 
described  in  [HSM83]  identify  necessary  and  sufficient  conditions  for  the  existence  of  a  C°° 
function  h(x)  such  that  r  =  n  and  provides  a  computational  approach.  The  necessary  and 
sufficient  conditions  for  (local)  EFL  are  nongeneric  and  not  likely  to  be  satisfied  in  general. 
In  the  sequel  we  show  how  that  essential  property  of  involutivity  of  the  /  and  g  vector 
fields  will  almost  never  be  satisfied  for  realistic  models  of  flexible  space  structures  due  to  the 
infinite  dimensional  nature  of  the  state  space. 

Byrnes  and  Isidori  [BI85]  describe  the  transformation  of  (2.3)-(2.4)  to  a  normal  form  in 
which  the  feasibility  of  PLF  control  can  be  assessed.  The  main  result  provides  the  existence 
of  a  diffeomorphic  transformation  of  coordinates  Th  :  3in  — ►  3?"  with  ( Th)(x )  (£,  z),  with 

the  state  partition  in  the  new  coordinates  £  £  9?n-r,  z  £  3fJr  and  inverse,  (2\)(*)  i— >  (£,  z),  so 
that  the  full  state  representation  in  the  new  coordinates  is 


where 


1  =  Fit,*), 

.  0  /r_l 

2  0  0 


Z  + 


0 

1 


H(£.*)  +  £(£,z)u], 


^,z)=  a(®)|x=fUi2) 


(2.20) 

(2.21) 


Definition:  The  zero  dynamics  of  the  input-output  model  (2.3)-(2.4)  are  given  by  the 
autonomous  system, 

(  =  F((,  0).  (2.22) 

And  the  system  is  locally  minimum  phase  if  the  the  diffeomorphic  transformation  Th  is 
defined  on  a  neighborhood  of  the  origin  and  (2.22)  is  asymptotically  stable  to  the  origin; 

£  =  o. 


2.3  PLF  Computations  for  Nonlinear  MIMO  Systems 
The  nonlinear  input/output  model  is 

x  =  f(x)  +  G(x)u  (2.23) 

y  =  h(x)  (2.24) 

where  x  £  9?n,  u,y  £  3f?m  and  f,gi  (resp.  yi)  for  i  =  1, . . .  ,m  are  smooth  vector  fields  defined 
on  3?"  (resp.  3?m).  For  notational  simplicity  we  write, 

G(x)  =  [g1(x),...,gm(x)]. 
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Nonlinear  PTant  Model 


Figure  2.2:  Partial  Linearizing  Feedback  via  Nonlinear  Inverse  Model 


Figure  2.3:  Partial  Feedback  Linearization  and  Zero  Dynamics 
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In  a  process  similar  to  the  previous  section  PLF  is  determined  by  transformation  of  the 
system  model  (2.23)-(2.24)  to  a  normal  form  in  which  we  identify  a  certain  ( mxm )  decoupling 
matrix  which  is  locally  nonsingular. 

The  process  begins  with  the  computation  of  an  appropriate  generalization  of  the  MIMO 
system  relative  degrees.  Let 

r-j  :=  min{fc  =  1,2, . . .  :  Lffi(L*-1(h,))  ^  0,  for  some  j  =  1, ...  ,m},  (2.25) 

the  itlx  characteristic  number  [FreT5j.  Each  r,-  is  then  the  minimal  relative  degree  of  the  set 
of  m  individual  output  responses  y{  obtained  from  each  input  Uj  for  j  =  1, . . .  ,m.  Let 

'  cti(x) 

A(x)  =  : 

**m(*) 

where  ai(x)  =  LJ(hi)  and 

‘  /3n(x)  •••  Pm(x) 

B(x)=  : 

flmm  (  *  ) 

where  fiij(x)  =  Lgj(Lr/~1(hi)).  Then  the  desired  normal  form  coordinates  are  z  €  9?r  where 
r  =  rj,  which  are  obtained  as 

Zi  > 

Z.2  (2-26) 
Zm  J 

with  each  z;  G  for  i  =  1, . . . ,  m  in  the  form, 


(2.27) 


Proposition:  Given  the  system  (2.23)-(2.24),  there  exists  a  diffeomorphic  transformation 
(T)(®)  =  (2,0  to  normal  form  coordinates, 

z  =  Az  +  E{A{z)  +  B(z)u)  (2.28) 

<  =  F(z,  f)  (2.29) 

where  A  =  diag{j4j, . . . ,  Am}  with 
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of  dimension  x  r;  for  each  t  =  1, . . . ,  m  and  E  is  an  r  x  m  matrix  with  elements  given  as, 


if  i  =  r,-  and  j  =  i 
otherwise 


Then  the  PLF  control , 

u=  -B~1(x){A(x)-v}, 

renders  the  v  y  input/output  response  in  linear  form, 


(2.30) 


i  =  Az  +  Ev,  (2.31) 

y  =  Cz.  (2.32) 


Definition:  The  system  (2.23)-(2.24)  (output  constrained)  zero  dynamics  are  given  by 

i  =  F(  0,0-  (2-33) 

Definition:  We  say  the  system  is  locally  minimum  phase  if  the  zero  dynamics  are  asymp¬ 
totically  stable  to  the  origin  (  =  0. 


Remark:  In  general  the  computation  of  the  normal  form  with  (2.29)  independent  of  u  is 
difficult  and  not  required  to  establish  the  minimum  phase  property.  Instead,  note  that  the 
system  zero  dynamics  are  just  the  dynamics  of  (2.23)-(2.24)  constrained  to  the  manifold 
Mh  C  3f?"  of  dimension  n  —  m  given  by, 

Mh  =  {iGr:  h(x)  =  o}. 

Proposition:  The  zero  dynamics  are  asymptotically  stable  if  and  only  if  the  system, 

*  =  /(*)  -  G(x)B~1(x)A(x),  *(0)  €  Mh 

is  asymptotically  stable  to  origin.  We  note  that  Mh  is  an  integral  manifold  lot  (2.23)-(2.24). 

2.4  Partial  Linearization  and  Variable  Structure  Control  Systems 

The  theory  of  Variable  Structure  (VS)  systems  addresses  the  design  of  control  laws  which  are 
discontinuous  functions  of  the  system  state.  VS  control  offers  practical  solutions  for  systems 
employing  actuators  which  can  be  efficiently  operated  in  bang-bang  and  other  discontinuous 
modes.  Our  interest  in  VS  control  for  rapid  slewing  of  multibody  systems  arises  from  the 
following  observations: 

1.  Design  methods  for  VS  control  for  output  regulation  have  been  shown  to  effect  an 
implicit  partial  feedback  linearization.  The  implicit  partial  linearization  is  achieved 
through  the  use  of  VS  control  requires  only  output  feedback. 
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2.  Direct  implementations  of  VS  control  laws  can  attain  a  level  of  robustness  to  plant 
model  assumptions  which  is  difficult  to  achieve  with  smooth  control.  Moreover,  ro¬ 
bustness  is  achieved  without  overly  conservative  restrictions  on  performance.  The 
limitations  of  robustness  for  VS  designs  are  clearly  related  to  the  minimum  phase 
conditions  for  PLF  design. 

3.  The  use  of  discontinuous  control  and  the  connections  between  VS  and  PLF  designs 
suggest  several  alternatives  for  integration  of  various  types  of  actuators,  including  both 
switched  and  continuous  modes  of  operation.  Integration  of  several  control  actuation 
systems  will  be  required  for  implementation  of  rapid  slewing  requirements  for  several 
candidate  large  space  platforms. 


The  general  theory  of  VS  control  design  is  well  known  and  we  will  not  attempt  to  present 
a  complete  description  of  its  scope.  For  details  see  the  survey  [DZM88].  However,  to  focus 
attention  on  the  concepts  we  seek  to  exploit  we  start  with  a  brief  description  of  the  basic 
ideas. 

VS  control  systems  utilize  high  speed  switching  control  to  drive  the  system  trajectories 
toward  a  specified  manifold  called  the  switching  surface.  Given  the  nonlinear  system  (2.23), 
the  VS  control  laws  are  of  discontinuous  type; 


Ui  = 


1 


«*(*). 


for  Si(x)  >  0 
for  Si(®)  <  0 


(234) 


with  s,(x)  =  0  smooth  switching  surfaces  chosen  in  the  state  space  for  each  i  =  1,...  ,m. 
The  design  approach  which  is  preferred  is  based  on  the  introduction  of  sliding  modes. 


Definition  (sliding  modes):  A  manifold,  At,,  consisting  of  the  intersection  of  p  <  m 
switching  surfaces,  s,(a:)  =  0,  with  the  property  that  a.ij  <  0  for  each  i  =  l,...,p  in 
the  neighborhood  of  almost  every  point  in  M,  is  called  a  sliding  manifold.  Under  these 
conditions  any  trajectory  of  the  system  which  enters  M,  remains  confined  to  the  manifold 
for  a  finite  length  of  time.  We  call  the  motion  on  M,  a  sliding  mode. 

VS  design  methods  involve  a  two  step  process:  1)  design  the  switching  surface  so  that  once 
sliding  is  achieved  the  natural  sliding  mode  achieves  design  objectives  such  as  regulation, 
stabilization,  etc.,  and  2)  design  of  discontinuous  control  laws  which  achieve  sliding  on  desired 
regions  of  the  switching  surfaces.  The  method  of  equivalent  control  is  a  popular  approach 
for  designing  the  switching  surface  to  achieve  desired  sliding  mode  dynamics. 

Given  the  system  (2.23)  and  a  manifold,  M,  =  {a:  €  :  s(x)  =  0},  with  s  :  5?"  — ► 

then  sliding  is  characterized  by  satisfaction  of  the  constraint  equations, 

s(x)  =  0,  i(sc)  =  0  (2.35) 

over  the  finite  time  interval,  <i  >  t  >  t2  where  s(i(ti))  =  0.  Note  that  a  sliding  mode  is  an 
instance  of  an  integral  manifold  for  the  closed  loop  system.  The  equivalent  control,  ue,,  is 
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the  control  required  to  maintain  the  system  trajectory  within  the  manifold  M.»  and  is  given 
by  the  condition, 

i  =  Vs(a:)  x  =  Vs(x){f(x)  +  G(x)utq}  =  0  (2.36) 

where  Vs(x)  =  ds/dx.  Under  the  assumption  that  det{Vs(a:)G(*)}  ^  0  for  x  £  M,  we 
have, 

utq  =  —  (Vs(a;)G(a:)]''1  Vs(x)/(®),  (2.37) 

and  the  motion  in  sliding  is  given  by, 

i  =  {I -G{x)[Vs{x)G(x)]-1Vs{x)}f(x),  s(x(t1))  =  0.  (2.38) 

Connections  between  the  design  of  VS  control  and  feedback  linearizing  control  have  re¬ 
ceived  considerable  attention  [FH87].  The  principal  focus  has  been  on  the  problem  of  synthe¬ 
sis  of  VS  designs  for  nonlinear  systems  of  the  form  (2.23)  using  (exact)  feedback  linearization. 

In  the  sequel  we  direct  attention  to  the  problem  of  output  regulation.  The  connection  we 
establish  with  PLF  design  also  illuminates  several  questions  relative  to  robustness  properties 
of  VS  designs  with  sliding  modes. 

Design  objective — Output  Regulation.  Given  the  system  (2.23)-(2.24)  where  y  indi¬ 
cates  a  set  of  regulated  outputs,  the  control  problem  is  to  drive  the  outputs  asymptotically 
to  zero. 

Since  output  regulation  problem  seeks  to  enforce  the  set  of  constraints 

hi(x)  =  0,  *  =  1, . . .  ,m, 

asymptotically  it  seems  reasonable  that  VS  design  could  be  employed.  However,  the  naive 
choice  s<(x)  =  h,(*)  leads  to  the  complication  that  in  general — in  fact,  most  often — [Vh(x)G(x)] 
is  singular  for  almost  all  a:  £  3?”.  The  approach  suggested  in  [KK89]  is  to  design  sliding  mode 
via  the  choice  of  switching  surfaces  relative  to  the  normal  form  coordinates  for  (2.23)-(2.24) 
as  given  by  (2.28)-(2.29). 

Proposition:  Given  the  system  (2.23)-(2.24)  obtain  the  diffeomorphic  transformation 
given  by  (2.26),  (2.27)  to  the  form  (2.28)— (2.29).  The  selection  of  switching  surface, 

s{z)  —  Kz  (2.39) 

with  K  an  m  x  m  constant  matrix,  solves  the  output  regulation  problem  if  sliding  can  be 
achieved.  In  sliding  the  equivalent  control  is, 

ueq  =  -B{x)~xKAz  -  B~\z)A{z)  (2.40) 

and  the  sliding  dynamics  are  given  by  the  r  linear  equations, 


z  =  [Ir-EK)Az,  Kz{  0)  =  0 
16 


(2.41) 
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Proof:  The  proof  is  given  in  detail  in  [KK89]  together  with  a  method  for  stabilization. 

By  defintion,  the  transformation  ( T)(x )  (,z,£)  is  invertible;  (T)(z,()  i->  x,  and  the 

switching  surface  s(z)  =  0  can  be  reflected  to  the  original  coordinates  s(®)  =  0.  The 
dynamics  in  the  x  coordinates  are  best  understood  in  terms  of  the  geometry.  Define  the 
m-dimensional  manifolds  M,  =  {x  £  :  s(x)  =  0}  and  Mh  =  {x  £  5?"  :  h(x)  =  0}.  The 

n  —  r  dimensional  manifold  Mz  =  {x  £  3?n  :  x  =  7(0, £)}  is  contained  in  both  Ma  and  Mh- 
Assume  that  in  some  neighborhood  T)  £  3 ?n  a  sliding  mode  exists  on  T)a  =  T>  Pi  M,  which 
is  assumed  nonempty.  Suppose  that  T>z  =  T>a  fl  Mz  is  nonempty.  Let  a  denote  a  bounded, 
stable  attractor  of  the  zero  dynamics  contained  in  Dz.  Assume  that  all  trajectories  in  T)z 
converge  to  a.  Then  if  the  initial  state  is  sufficiently  close  to  V,  the  trajectory  will  eventually 
reach  T>,  and  sliding  will  occur.  Clearly  the  stability  of  the  attractor  is  critical  to  the  stability 
of  the  overall  design.  In  the  sequel,  we  establish  conditions  for  output  regulation  of  multibody 
systems  which  guarantee  that  the  zero  dynamics  have  well  defined  local  equilibria  so  that 
linear  stability  analysis  of  the  zero  dynamics  is  appropriate.  We  remark  that  the  problem  of 
establishing  estimates  for  the  domain  of  attraction  in  the  zero  dynamics  is  an  open  question. 

2.5  Robust  Stabilization  of  Nonlinear  Systems 

For  practical  implementation  of  PLF  compensation  various  researchers  have  focused  atten¬ 
tion  on  conditions  which  guarantee  robust  stabilization  of  the  nonlinear  system  (2.3)-(2.4) 
with  feedback  transformation  of  the  form  (2.11)  by  introduction  of  linear  feedback  v  =  Kz. 
A  brief  survey  of  the  wide  range  of  methods  which  have  been  proposed  is  given  in  [BBKA88]. 
In  this  section  we  focus  attention  on  a  ubiquitous  assumption  in  most  of  the  work  on  robust 
stabilization  of  nonlinear  systems  based  on  PLF  compensation. 

Consider  the  usual  case  for  engineering  design  where  the  open  loop  system  dynamics  for 
(2.23)-(2.24)  is  given  by  a  nominal  model  of  the  form 

x  —  f°(x)  -f  G°{x)u  (2-42) 

V  =  h°(x),  (2.43) 

where  f° , g°  are  C°°  vector  fields  for  i  =  l,...,m,  defined  on  a  manifold  M  £  9?n,  with 
f°{ 0)  =  0.  We  assume  the  (true)  system  response  can  be  modeled  via  a  perturbation  of  the 
vector  fields; 

/  =  /°  +  A/,  G  =  G°  +  AG 

with  A/,  Agi  each  C°°  defined  on  M  and  A/(0)  =  0.  In  [AB88]  a  detailed  analysis  is  given 
leading  to  sufficient  conditions  on  A/,  AG  which — together  with  the  assumption  that  (2.42) 
is  feedback  linearizable — guarantees  that  (2.23)  is  also  linearizable.  The  conditions  given  are 
less  restrictive  then  the  usual  structure  matching  conditions  [AB87].  The  structure  matching 
conditions  also  play  a  role  in  establishing  conditions  for  robust  stabilization  and  we  repeat 
them  for  convenience. 
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The  structure  matching  conditions.  Under  the  assumption  that  the  perturbation  vec¬ 
tor  fields  satisfy, 

A/,  Ag{  6  A  (2.44) 

where  A  =  Sp{<7i, . . . , pm }  then  any  such  model  (2.23)  is  exactly  linearizable — in  fact,  by  the 
same  diffeomorphic  transformations.  The  above  conditions  are  equivalent  to  the  statement 
that  the  perturbations  to  the  vector  fields  can  be  factored  as: 

A /(*)  =  G°(x)d,(x)t  (2.45) 

A  G(x)  =  G°(x)Dg{x).  (2.46) 

The  importance  of  the  structure  matching  conditions  in  establishing  robust  stability  is 
that  under  these  conditions  the  model  uncertainty — after  application  of  PLF  compensation- 
can  be  equivalently  represented  by  a  perturbation  (or  disturbance)  at  the  compensated  sys¬ 
tem  inputs.  This  facilitates  the  design  of  disturbance  rejection  techniques  using  either  ex¬ 
plicitly  nonlinear  control  designs  such  as  Gutman  [GL76]  or  linear  control  design  such  as  in 
[Kra87].  To  see  this  we  summarize  the  construction  under  the  assumption  that  h(x)  =  h°(x). 
Substitute  the  PLF  compensation  (2.30)  obtained  for  the  nominal  model  (2.42)-(2.43)  into 
the  model  for  the  true  system  (2.23)  to  obtain, 

i  =  f°  +  Af  +  {G°  +  AG)B~l{v  -  A) 

=  f°  +  G°Ad,  +  G°(I  +  Dg)B~l[v  -  A] 

=  f° +  G°B-1[v-A}  +  G°[df  +  DgB{v  -  A)}.  (2.47) 

The  model  after  nominal  PLF  compensation  can  (by  the  above  assumptions)  be  transformed 
to  z-coordinates  defined  in  (2.26)-(2.27)  to  obtain, 

z  =  Az  +  E{v  -f  7]} 

i  =  F(*,a 

where 

<1  <*,{,»)  =  (2.50) 

=  [</W  -  |>=r_,(i)  +  [£>,(*)0(*)-v]  |i=r_1(>) .  (2.51) 

Thus  it  is  clear  that  robust  stabilization  must  address  the  disturbance  rejection  of  the  class 
of  input  disturbances  dv(t)  which  bound  the  model  error;  ||j/||  <  {|<fv ((.  An  effective  design 
approach,  in  the  case  when  the  nominal  model  is  exactly  feedback  linearizable,  is  given  by 
Spong  and  Vidyasagar  [SV87].  There  approach  utilizes  an  Loo  stabilization  criterion  and 
obtains  a  linear,  time-invariant  feedback  control  for  the  v-input. 


(2.48) 

(2.49) 
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Remarks?  The  design  methods  for  robust  stabilization  of  nonlinear  systems  available  in  the 
literature  are  almost  exclusively  based  on  the  structure  matching  conditions.  By  reflecting 
the  model  uncertainty  to  the  system  inputs  (after  nonlinear  feedback  compensation)  they 
can  employ  either: 

1.  linear  compensator  design  which  seeks  to  reduce  loop  gains  consistent  with  bounded 
model  uncertainty,  or 

2.  nonlinear  switching  mode  compensator  design  which  seeks  to  over-bound  input  distur¬ 
bances  by  high  gain  implementations  using  fast  switching  control. 

Both  methods  result  in  essentially  conservative  designs  since  the  worst  case  bounds  on  the 
input  disturbances  must  be  assumed. 

In  [Akh89]  the  basis  for  robust  stabilization  of  nonlinear  systems  is  considered  further 
and  new  results  are  obtained  for  the  case  of  parametric  model  uncertainty.  The  new  con¬ 
trol  laws  obtained  in  [Akh89]  employ  basic  constructions  of  adaptive  control  in  the  context 
of  feedback  linearization.  The  results  show  that  robust  stabilization  can  be  obtained  under 
much  less  restrictive  conditions  than  the  structure  matching  conditions.  Significantly,  feed¬ 
back  linearization  plays  the  role  of  enforcing  linearity  for  subsequent  control  loop  designs. 
To  the  extent  that  this  can  be  achieved  in  practical  applications  it  can  enhance  reliability 
and  repeatability  thus  achieving  improved  performance  prediction— an  important  feature 
for  space-based  systems.  Standard  constructs  in  adaptive  control  can  then  be  applied  to 
enhance  the  robustness  of  feedback  linearization  with  model  uncertainty.  We  emphasize 
that  the  constructions  described  below  are  new  and  offer  stability  results  for  the  nonlinear 
system  under  very  general  assumptions.  In  the  next  few  paragraphs  we  briefly  review  some 
significant  aspects  of  these  results. 

Robust  Stabilization  of  Nonlinear  Systems  by  Adaptive  Methods.  Again,  starting 
with  the  system  model  in  the  form  (2.23)  we  assume  the  model  uncertainty  can  be  represented 
by  parametric  dependence  of  the  vector  fields  so  that  the  model  has  the  form, 

x  =  f{x,0)  +  G(x,9)u,  (2.52) 

with  9  a  p-vector  of  unknown  parameters.  We  assume  that  for  every  9  £  Be,  a  closed, 
compact  neighborhood  of  the  nominal  parameter  90,  f  and  gi,  for  i  =  l,...,m,  are. C°° 
vector  fields  and  f(Q,9)  =  0  for  9  £  Be .  The  nominal  design  model  is  characterized  by  the 
set  of  nominal  parameters  and  we  take  f(x)  =  f(x,9a),  G(x)  =  G(x,9a). 

The  following  assumptions  are  used  in  [Akh89]  to  establish  a  robust  stabilizing  controller 
for  the  nonlinear  system. 

Assumption  1:  The  nominal  system  is  exactly  feedback  linearizable.1 

*The  extension  r  f  the  results  described  below  to  the  case  of  stabilization  by  PLF  is  in  progress. 
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It  is  important  to  note  that  assumption  1  is  applied  only  at  the  nominal  plant  model 
with  fixed  and  known  parameters  90.  This  is  in  contrast  to  the  structure  matching  conditions 
which  are  almost  universally  assumed  in  the  current  literature. 

Assumption  2:  For  any  x  €  U  C  5?"  and  9  (E  Bo, 

&gi{x,9)  e  Sp  {gi{x,0„)t. . .  ,gm(x,90)} 

for  z  =  1, . . . ,  m.  This  assumption  implies  that  there  exists  an  m  x  m  matrix  valued  function, 
D(x,9),  with  smooth  elements  such  that  A G  =  G°D. 

Assumption  3:  Either  there  exists  an  m  x  m,  strictly  positive  definite  matrix,  K  such 
that, 

Vx  eU,V9e  Bo,  0  <  D{x,9)  <  K, 
or  there  exist  a  K,  negative  definite,  such  that, 

Vx  G  V,  \/9  e  Bo,  K  <  D{x,  9)  <  0. 

This  assumption  (in  various  forms)  is  typical  in  adaptive  control  stability  analysis  and 
design.  It  says  that  the  “sign”  of  this  term  must  be  definite  and  known  a  priori. 

The  design  approach  is  natural  and  begins  with  the  transformation 

2  =  T(x,0o), 

to  normal  form  and  choice  of  feedback  linearizing  compensation  (2.30).  The  v  control  is 
chosen  in  two  parts.  First,  for  the  nominal  design  and  performance  objectives  we  find 
v  =  Fz  where  A  -f  EF  in  (2.31)  is  a  stable  matrix.  In  this  case,  there  exists  a  unique, 
positive  definite,  symmetric  solution  to  the  Lyapunov  equation, 

{A  +  EF)T P  +  P(A  -t  EF)  =  -I. 


The  design  for  the  nominal  model  is  now  modified  by  the  introduction  of  adaptation. 
The  control  law  obtained  is  described  by  the  following  equations; 


u  =  —B  1(®){A(x)  —  v), 

(2.53) 

2  =  T(x,90), 

(2.54) 

v  =  Fz  +  C(z,0)9, 

(2.55) 

9  =  -T-1Ct{z,9)EtPz, 

.  (^  56) 

C(z,9)  =  -pB{z)GT(z,90)Pz9T. 

(2.57) 

In  these  control  laws  the  pxp  matrix  T  >  0  is  thu  “adaptation  gain” 

which  is  chosen  so  that 

the  Lyapunov  function, 

V  =  zTPz  +  6tT9 

is  positive  definite  for  all  x  €  U  and  p  is  a  positive  scalar.  The  result  established  in  [Akh89] 
is  the  following. 
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Theorem:  The  adaptive  control  with  feedback  linearizing  transformation  is  asymptotically 
stable  to  the  origin  x  =  0  and  9-^9,  asymptotically,  if  for  x  6  U  and  9  £  B$  there  exists 
cj  >  1  and 

2||  TT(x)P9(x,e )  -  nTDy  ||»||!  ||  <  c2||r(x)||J, 

where 

9(x,  9)  =  [A f(z,9)  +  A G{z,  0){A(z)  +  B~lFz}\ 

and 

y  -  Gt(z,0o)Ptz. 


Proof:  [Akh89]. 
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y,o 


x,<f> 


Figure  3.1:  Standard  Coordinate  Frame  for  Modeling 

3  Nonlinear  Dynamics  of  Multibody  Systems  with  Flexible  In¬ 
teractions 

In  this  section  we  describe  the  basis  for  the  formal  development  of  a  class  of  evolution  mod¬ 
els  for  multibody  systems  with  elastic  interactions.  We  follow  the  approach  suggested  by 
Baillieul  and  Levi  [BL87].  As  will  be  seen,  the  formulation  captures  the  essential  evolution 
dynamical  structure  of  the  system  without  requiring  detailed  knowledge  of  its  internal  con¬ 
figuration.  As  such  this  framework  provides  a  consistent  modeling  approach  for  developing 
a  hierarchy  of  models  with  increasing  internal  complexity  and  fine  structure.  The  main  idea 
is  to  isolate  a  “primary  body”  and  to  attach  a  reference  frame  to  it  at  a  convenient  loca¬ 
tion  for  measuring  attitude  and  displacement  dynamics.  The  motion  of  all  other  spacecraft 
components  will  then  be  measured  relative  to  this  frame. 

Throughout  this  report  we  will  use  the  notational  conventions  given  in  Table  3.1.  Con¬ 
sider  a  reference  frame  fixed  in  the  primary  body,  with  origin  located  by  the  position  vector 
R  €  S?3  and  angular  orientation  denoted  by  L  G  SO( 3),  both  relative  to  a  fixed  inertial  frame 
(see  Figure  3.1).  L  can  be  parameterized  by  the  Euler  angles2  ip,8,<p  representing  sequential 
rotations  about  the  axes  3,2,1,  respectively: 

cos  0  cos  ip  cos  8  sin  6 

L  =  sin  <p  sin  0  cos  ip  —  cos  <f>  sin  ip  sin  <p  sin  8  sin  ip  4-  cos  <p  cos  ip 
cos  <p  sin  0  cos  ip  4-  sin  <p  sin  ip  cos  <p  sin  8  sin  ip  —  sin  <p  cos  ip 

2  We  use  the  so-called  NASA  standard  or  321  convention  [Gol82]. 


—  sin  0 

sin  <p  cos  0  (3-1) 

cos  <p  cos  8  _ 
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Notation 


Xj,  t  —  1,2,... 


x  = 


dx 

dt 


Xz(t,z)  =  fz  (z,t) 


6x 

x  x  y 


Explanation 

element  of  a  vector  x 

transpose  of  vector  x 

time  differentiation 

partial  differentiation 

natural  (Hilbert  space)  inner  product 

natural  (Hilbert  space)  norm 

differential  variation 

vector  cross  product  of  x  and  y 


Table  3.1:  Standard  Notation  for  Lagrangian  Mechanics 


A  fundamental  kinematic  relationship  is 

L(t)  =  -fi  (t)L(t) 


(3.2) 


where 


‘  0 

—(*>3 

U>2 

n  = 

U>3 

0 

-wi 

.  —(^2 

U>1 

0 

(33) 


and  u)  =  is  the  primary  body  (inertial)  angular  velocity  as  measured  in  the 

body  coordinates. 

Define  (  £  §?3  as  £  =  (ip,  9,  (p)T .  Then  an  equivalent  relation3  is 


*  =  r(£K 


1  0  —  sin  9 

0  cos  ip  cos  6  sin  ip 

0  —  sin  ip  cos  6  cos  ip 


(34) 


The  body  frame  position  and  orientation  can  be  characterized  as  a  point  in  the  special 
Euclidian  group,  SE( 3,3?),  each  element  of  which  can  be  represented  by  a  matrix 


X  = 


'  LT  R 
0  1 


(3.5) 


The  positions  of  all  other  elements  of  the  system  are  measured  relative  to  the  primary  Body 
frame.  We  identify  each  particle  (or  element),  P,  by  its  “undeformed”  position,  z,  in  the 
primary  body  frame.  Let  u(z,  t)  denote  the  deformed  position  of  P.  Furthermore,  we  fix 
a  coordinate  system  in  each  particle  with  origin  at  u(z,t)  and  aligned — in  the  undeformed 
state — with  the  body  axis  coordinates.  Let  l(z,  t)  £  50(3)  denote  the  orientation  of  P  in  the 

3Equation  (3.4)  is  essential  to  the  analytic  framework  for  multibody  modeling.  Alternate  parametrisations 
of  50(3) — such  as  the  Cayley-Rodrigues  parameters  [Dwy84] — can  be  used  to  advantage  and  the  general 
form  remains  intact. 
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deformed  state  as  measured  in  the  primary  body  coordinates.  Note  that  in  the  undeformed 
state 


(unde  formed  — ' 


1  0  0 


0  0  1 


and  for  small  relative  motions4 


■  1  ^  —$ 

(small  =  1  <f>  .  (3-7) 

0  -<f>  1 


The  inertial  coordinates  U(z,t )  of  a  particle  P  can  be  obtained  from  the  body  coordinates 
u (z,  t)  via  the  relation 

U(z,t)  =  Ltu{z,  t)  +  R.  (3.8) 

Note  also  that 

_  f  LT{t)m  m 
X{t)  -  \  0  0 

Also,  a  direct  computation  yields 

—U  =  LT(t)[Uu  +  u]  +  R{t)  (3.10) 

at 

The  kinetic  energy  of  the  system  can  be  written  in  terms  of  the  generalized  coordinates 
q  =  (£,  R,  u)  in  the  form 

T(q,q)  =  iJs\\U\\2dm  =  i  JjL^Uu  +  ut]  +  R\\2dm 

=  \  Js\\[nu  +  ut]\\2  +  2{nu  +  ut,LR)  +  ||J2||**n  (3.11) 

where  S  denotes  that  the  integral  is  to  be  taken  over  the  entire  system. 


3.1  Lagrange's  Equations  for  Continuum  Dynamics 

The  formalism  of  Lagrangian  dynamics  begins  with  the  identification  of  the  configuration 
space,  i.e.  the  generalized  coordinates,  associated  with  the  dynamical  system  of  interest. 
Once  the  configuration  manifold,  M,  is  specified  we  have  the  natural  definition  of  velocity 
at  a  point  q  G  M  as  an  vector,  q,  in  the  tangent  space  to  M  at  g,  often  denoted  TqM.  We 
then  define  the  state  space  as  the  union  of  all  points  g  6  M  along  with  their  tangent  spaces, 
the  so-called  tangent  bundle  (c.f.  [AM78,  Arn78])  TqM.  The  evolution  of  the  system  in 
the  state  space  is  characterized  using  Hamilton’s  principle  of  least  action  by  the  definition 

4The  assumption  of  small  relative  motions  is  useful  in  the  continuum  framework  to  model  large  displace¬ 
ments.  This  idea  suggests  several  alternative  approximations  useful  for  modeling  structural  flexure  dynamics 
and  is  discussed  in  detail  in  $  3.3. 
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of  a  Lagrangian  L(q,q)  :  M  x  TqM  — ►  $?.  Hamilton’s  principle  says  that  the  motion  of  a 
dynamical  system  between  times  t\  and  t%  is  a  “natural”  motion  if  and  only  if 


r1 1 

6  L  dt  =  0, 
Jti 


(3.12) 


or — accounting  for  the  presence  of  external  generalized  forces,  Q — in  its  generalized  form; 

(3.13) 


+  QTSq)  dt  =  0. 
Jt, 


For  distributed  parameter  systems  (DPS),  special  care  is  required  to  properly  characterize 
the  configuration  space  for  modeling  the  system  motions.  The  principal  reasons  for  this 
fact  follow  from  the  application  of  the  models  obtained;  viz.,  the  study  of  time  evolutions 
subject  to  control  forces.  First,  control  systems  will  inevitably  involve  the  implementation 
of  feedback  and  we  are  therefore  immediately  concerned  with  stability.  An  appropriate 
notion  of  stability  is  central  to  the  design  of  feedback  control  systems.  For  Lagrangian 
systems  the  natural  definition  of  stability  is  implicit  in  the  structure  of  the  state  space 
which  for  DPS  is  a  function  space  and  care  must  be  exercised  that  the  construction  (and 
assumptions)  of  the  state  space  are  consistent  with  the  engineering  control  problem.  Second, 
it  is  often  necessary  to  define  finite  dimensional  approximations  to  DPS  for  a  variety  of 
reasons  including  computer  simulation.  Again,  our  primary  concern  is  in  approximating  the 
time  evolution  under  the  influence  of  control.  As  we  will  make  clear  in  the  following  section, 
the  formulation  of  such  models  in  a  consistent  way  is  inherently  bound  to  the  definition  of 
the  configuration  space.  In  this  section  we  confine  our  discussion  to  the  configuration  space 
for  continuous  systems  with  one  spatial  dimension. 

The  generalized  coordinates  are  chosen  so  that  all  “nonworking”  or  geometric  constraints 
on  the  motion  are  eliminated.  This  is  the  key  to  the  utility  of  the  Lagrange  formalism 
for  constructing  the  equations  of  motion.  In  the  case  of  DPS  any  “geometric”  boundary 
conditions  (which  we  will  denote  Q)  are  therefore  included  as  part  of  the  definition  of  the 
configuration  space.  All  other  boundary  conditions  necessary  to  complete  the  Euler- Lagrange 
equations  result  from  the  application  of  Hamilton’s  Principle,  (3.12)  or  (3.13).  These  are  the 
“natural”  boundary  conditions  (denoted  M). 

An  essential  part  of  the  definition  of  the  configuration  space  in  the  infinite  dimensional 
case  is  the  specification  of  the  norm.  Although  all  norms  are  equivalent  in  finite  dimensions, 
this  is  certainly  not  the  case  in  infinite  dimensions.  We  briefly  summarize  the  main  issues. 
Consider  functions  v(z)  defined  on  the  domain  z  €  [0, 1]  and  let  Drv(z)  denote  the  rth 
derivative  with  respect  to  z.  We  denote  by  Hp  the  completion  of  the  set  of  the  set  of 
functions  with  p  continuous  derivatives  and  which  satisfy 

IMG  -  /o'{|£”M*)P  +  ■■  •  +  IvMI’H*  <  oo  (3.14) 

These  are  the  Sobolev  spaces  [Lio71] .  Equivalently,  Hp  consists  of  those  functions  whose  first 
p  derivatives  belong  to  the  Hilbert  space  of  square  integrable  functions.  Note  that  x  €  H* 
implies  *  e  H'~l  for  i  =  1, 2,  — 
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Let  Hg  denote  the  completion  of  the  set  of  functions  satisfying  (3.14)  as  well  as  a  pre¬ 
scribed  set  of  boundary  conditions  designated  Q.  It  is  not  necessarily  true  that  all  of  the 
functions  in  this  new  space  satisfy  the  boundary  conditions.  The  reason  for  this  is  that  an 
arbitrary  sequence  of  functions,  all  satisfying  the  given  boundary  conditions,  may  converge 
to  a  function  which  does  not  satisfy  the  boundary  conditions.  However,  the  following  propo¬ 
sition  is  true.  Suppose  the  boundary  conditions  Q  involve  derivatives  of  order  s  and  none 
higher.  Then  all  of  the  functions  in  Hg  satisfy  the  boundary  conditions  provided  p  >  s. 
Thus,  a  consistent  definition  of  the  configuration  space  is  obtained  if  the  specified  norm  is 
compatible  with  the  geometric  boundary  conditions. 

Hamilton’s  principle  may  be  used  to  derive  the  Euler-Lagrange  equations  and  the  natural 
boundary  conditions.  The  Euler-Lagrange  equations  are  to  be  solved  along  with  boundary 
conditions  B  =  Q  U  N.  In  general,  the  Lagrangian  will  involve  derivatives  with  respect  to  z 
of  order  p  and  the  Euler-Lagrange  equations  will  involve  derivatives  of  order  2 p.  In  finding 
solutions  q(t)  we  seek  “weak”  (sometimes  called  generalized  or  distributional5)  solutions  in 
Hg  which  satisfy  Hamilton’s  principle  or  “strong”  (pointwise,  genuine  or  classical)  solutions 
in  Hgp  which  satisfy  the  Euler-Lagrange  equations.  The  results  are  equivalent  (in  Hg)  when 
both  problems  have  solutions.  The  Euler-Lagrange  equations  may  be  given  the  interpretation 
of  an  evolution  equation  as  we  sill  describe  below. 

Finite  Dimensional  Approximation  and  Computer  Simulation  Finite  dimensional 
approximations  to  the  system  dynamics  may  be  obtained  by  seeking  an  approximate  solu¬ 
tion  to  the  Euler-Lagrange  equations  or  to  Hamilton’s  principle  directly.  The  latter  has  the 
advantage  that  solutions  are  to  be  sought  in  a  larger  space  of  admissible  functions  which  pro¬ 
vides  a  wider  choice  of  approximating  functions.  Perhaps  unexpectedly,  this  turns  out  to  be 
of  fundamental  significance  in  developing  numerical  solutions  to  the  required  evolution  dy¬ 
namics  and  for  computer  simulation.  Furthermore,  important  links  to  the  system  physics  are 
retained  through  this  modeling  process.  These  observations  appear  consistent  with  many 
standard  engineering  methods  which  introduce  approximations  to  continuous,  distributed 
system  dynamics  by  discretization  of  the  variational  problem  underlying  the  Lagrangian 
dyanmics  [Mei67].  Indeed,  this  is  the  basis  for  the  Finite  Element  Method  (FEM)  for  evolu¬ 
tion  dynamics  described  in  [SF73].  The  simulation  models  developed  in  this  study  are  based 
on  finite  dimensional  approximation  using  collocation  by  splines  [Aga84,  Sta79,  Pre75].  Fur¬ 
ther  details  of  the  method  will  be  given  in  a  later  section  with  examples. 

Next  we  consider  some  simple  continuous  systems  arising  in  structural  mechanics  which 
will  illustrate  the  evolution  modeling  setup  described  above. 

5There  are  several  approaches  to  developing  the  notion  weak  and  strong  solutions  and  many  good  reasons 
for  doing  so.  Thus,  the  proliferation  of  terminology  carries  with  it  sometimes  subtle  distinctions  (c.f.  [SF73, 
RM57,  Sta79j). 
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Figure  3.2:  Simple  Cantilevered  Beam 
3.1.1  Example:  Simple  Cantilevered  Beam 

Consider  the  cantilevered  beam  undergoing  small  transverse  motions  confined  to  the  plane. 
The  beam  is  excited  by  a  concentrated  force,  /(£),  and  moment,  m(£),  applied  at  the  point 
z  =  a  €  (0, 1). 

Timoshenko  Model.  Each  cross  section  undergoes  a  displacement  t](z,  t)  and  a  rotation 
<f>(z,  t).  These  are  the  generalized  coordinates.  The  geometric  boundary  conditions  are 

G  :  77(0, 1)  =  0  and  <£( 0, 1)  =  0. 

Thus,  the  appropriate  configuration  space  is  Hg.  The  Lagrangian  is 


and  the  virtual  work  8W  =  QT6q  due  to  the  external  forces  is 

8W  =  f  {f(t)8(z  —  a)8r\  +  m(t)8(z  —  a)8<f>}dz, 

Jo 

where  p  is  the  mass  density,  A ,  the  cross  section  area,  /,  the  moment  of  inertia,  E,  the 
modulus  of  elasticity,  and  kG,  the  effective  shear  modulus.  Upon  application  of  Hamilton’s 
principle,  we  obtain  the  partial  differential  equations 

and  the  natural  boundary  conditions 

-“pt2-*)- - 
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Thus,  we  have  the  evolution  equation 


pA  0 
0  pi 


1(5)- 


d_  f  kGA§; 


dz  kGA  ei 


■kGA  (  tj  \  _  (  f(t)  \ 
j;kGA  V  <f>  )  V  m(t)  J 


6(z  —  a) 


where  we  interpret  [q(-,t),<l>(-,t)]T  as  an  element  in  HI 


Bernoulli-Euler  Model.  Suppose  that  we  consider  the  same  situation  with  the  additional 
Bernoulli- Euler  assumptions  [BK89].  These  are 

1.  rotational  inertia  is  negligible,  pi  — >  0, 

2.  shear  deformation  is  negligible,  i]z  —  <f>  — >  0. 

The  deformed  beam  configuration  is  completely  specified  by  t ](z,  t).  The  geometric  boundary 
conditions  are 

Q  :  ?/(0, 1)  =  0  and  1)  __  q 

oz 

Notice  that  the  appropriate  configuration  space  is  Hg .  The  simplified  Lagrangian  is 


and  the  virtual  work  expression  also  simplifies  to 


=  j  ^f{t)8(z  -  a)8r)  +  m(t)6(z- a)8  jj 
=  J  |/(t)5(z  —  a)  +  m(t)8~1(z  —  a)}  8rj  dz. 


(3.15) 


The  evolution  equation  is 


=  f(t)&{z- a)  +  m(t)6  1{z -a), 


which  is  to  be  interpreted  on  Hg  with 


dz9  dz 2 
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3.2  Lagrange’s  Equations  for  Multibody  Systems  with  Elastic  Interactions 

In  addition  to  the  kinetic  energy,  T(q,q),  we  assume  that  a  potential  energy  function  V(q ) 
is  also  available.  Then  Lagrange’s  equations  take  the  form 

d_  dT  _  dT  dV_ 
dt  di  ~  Q(' 


±  dT  _dT  dV 

dt  dR  ~  dR  +  dR  ~  Qr' 
&_  ST  _  ST  8V 

dt  8u  8u  ^  8u 


(3.16) 


where  the  generalized  forces  are  defined  in  terms  of  the  virtual  work  expression ; 

8W  =  Qjd(  +  QldR  +  QtJu.  (3.17) 

Now,  we  define  the  system  angular  momentum  with  respect  to  the  origin  of  the  body  frame 

[(a;  x  u  +  ut)  +  LR]dm  =  a  -f  J  u  x  LR  dm.  (3.18) 

With  some  calculation6  these  equations  reduce  to 

System  Angular  Momentum 

rT(|)[d  -fa;  x  a]  +  J^u(z,t)  x  LR  dm  =  (3.19) 

System  Linear  Momentum 


Is^dtLT(Ut  +  w  x  w)  +  dm  =  Qr~ 

8V 

+  w  x  (w  x  «)  +  «  x  tt  +  2w  x  uj  +  LR  =  Qu  — — . 

8u 


(3.20) 

(3.21) 


Equivalently  we  obtain 

rr(£)[/w+o> X  Iw]+mcxLR+TT(t)D  f  ux  [wxu+ut]dm+  f  ux  LRdm  =  Q^  —  (3.22) 

J  S  J  S 


mR  -f  mZrr[a;  x(wxc)  +  wxc]  +  LT[D2u  +  LR\dm  =  Qr 


D2 u  -f  LR  —  Qu  — 


(3.23) 

(3.24) 


where  the  operator  D  is  defined  by 


£(•):=  j;(') +  “'*(•) 

Note  that,  in  applications,  the  integrals  in  (3.19)-(3.21)  or  (3.22)-(3.21)  would  not  be 
evaluated  directly.  Instead,  they  are  to  be  replaced  by  momentum  expressions  in  terms  of 
an  appropriate  choice  of  generalized  coordinates. 


8 Appendix  A  contains  some  identities  useful  for  these  calculations. 
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3.3  Geometrically  Consistent  Reductions  and  Modeling  of  Flexible  Structures 

Recently  attention  has  been  drawn  to  gross  qualitative  properties  of  certain  dynamic  sim¬ 
ulation  models  for  multibody  aerospace  systems  with  structural  flexure  which  have  been 
obtained  by  “conventional  methods”.  Such  modeling  methods,  which  often  employ  approx¬ 
imations  involving  linearization  and  truncation  of  position  and  velocity  expressions  refer¬ 
enced  to  body  fixed  frames,  can  lead  to  gross  qualitative  errors  in  predicted  dynamic  re¬ 
sponses.  Examples  include  inconsistent  or  incorrect  predictions  of  dynamic  stiffening  or 
softening  when  the  body  frame  is  subject  to  arbitrary  large  motion  in  inertial  space.  In 
certain  examples  dynamic  predictions  suggest  unstable  responses  inconsistent  with  physi¬ 
cal  experiments.  Such  errors  in  qualitative  behavior  of  models  can  be  best  understood 
by  a  careful  analysis  of  the  geometry  of  the  deformations  allowed  under  certain  approxima¬ 
tions  [SVQ86,  VQS87,  SVQ88,  Pos88].  Several  computer  codes  for  multibody  modeling  simu¬ 
lation  of  Control-Structure  Interaction  (CSI)  effects  such  as  DISCOS,  TREETOPS,  NB0D2, 
ALLFLEX,  are  widely  available  and  used  in  the  aerospace  industry.  The  approximations 
used  implicitly  in  these  codes  have  been  repeatedly  called  in  to  question  by  representatives 
of  both  industry  and  academia  [KRB87,  BD90]. 

In  this  section  we  briefly  describe  the  geometry  of  large  deformations  from  the  viewpoint 
of  continuum  Lagrangian  methods.  We  outline  some  methods  of  approximation  and  reduc¬ 
tion  of  models  which  are  “geometrically  consistent”.  Two  issues  arise  with  respect  to  the 
geometry  of  large  deformations,  the  first  is  the  notion  of  strain  and  its  definition  in  terms 
of  local  material  deformations.  The  second  arises  when  structures  exhibit  large,  transverse 
global  deformations  but  one  chooses  to  make  certain  kinematic  assumptions  which  restrain 
axial  deformation. 

3.3.1  Definition  of  Strain 

Consider  a  continuous  medium  which  undergoes  a  deformation.  Let  I  6  S3  denote  the 
coordinates  of  a  particle  P  lying  on  a  curve  Co  in  the  undeformed  state.  Suppose  the 
curve  C0  deforms  to  the  curve  C.  After  deformation  the  location  of  the  particle  P  on  C  is 
denoted  *  €  5n  .  If  we  consider  the  particles  to  be  infinitesimal  rathe  than  points,  then  the 
characteristic  length  of  the  curve  Co  (resp.  C)  is  the  arc  length  ds0  (resp.  da).  Similarly, 
the  respective  infinitesimal  area  and  volume  elements  are  proportional  to  da2  and  da3.  We 
may  write 

dal  —  dX*  +  dX l  -I-  dX\  and  da 2  =  dx\  -f  dx\  +  dx\. 

Taking  the  viewpoint  of  Lagranian  mechanics  [Sok56,  Fre66]  we  identify  the  location  of  each 
point  P  relative  to  its  nominal  location  in  the  undeformed  state  of  the  structure;  i.e.,  the 
coordinates  are  functionally  dependent;  *  =  x(X).  Then  the  differential  volume  element  can 
be  given  as, 

-.in— {[unm}-  <3-2s> 
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We  would  like  to  express  strain  as  a  measure  of  deformation  per  unit  area.  We  thus  define 
strain  in  terms  of  the  area  deformation, 

ds2  -  ds2  :=  2  dXt  E  dX. 

where  E  is  the  (Lagrangian)  strain  tensor ,  which  satisfies, 

„  „  f  dx  1 4  f  dx  1  ... 


2  E  = 


dX\  \dX 


(3.26) 


It  is  convenient  to  express  the  strain  in  terms  of  the  relative  displacement  by  introducing 
u(X)  which  satisfies  the  relation; 


so  that 


IE  = 


x{X)  =  u(X )  +  X, 

du  1  *  [  du  1 

+  -»3  7TT7  +  -*3  _  1 3- 


From  (3.27)  we  obtain  the  components  of  the  strain  tensor,  tij  for  i,j,  k  =  1,2,3  as, 


(3.27) 


(3.28) 


and  for  i  ^  j  ^  k, 


n  __  9ui_  dui  fdv±  dui_  du^  dvj_  dy*  du*) 

lJ  dXj  dXi  \  dXi  dXj  dXi  dXj  8Xi  dXj  j '  [  ) 

When  deformations  are  small — i.e.,  when  u(X)  and  its  partial  derivatives  are  negligible — 
the  second  order  terms  in  the  strain  tensor  can  be  neglected  and  we  recover  the  infinitesimal 
strain  definitions  [Sok56,  Fre66], 

dUi  / 


dtt;  dui 

*«  =  WSM'  hti^ 


(3.30) 

(3.31) 


Remarks  on  Extension:  Consider  a  line  segment  Co  which  lies  in  the  Xi  axis.  At  any 
point  P  in  Co  an  element  length  is  dso  =  dX\.  Thus,  using  the  definition  of  strain  we  have, 

ds2 -ds2  =  2dXt  EdX  =  2endX2. 

The  extension  of  this  element  under  deformation  is  given  by, 

E  =  (ds  —  dso)/ds0  =>  ds  =  (1  -j-  E)dso. 

Then 

ds2  -  ds\  =  [(1  -f  E)2  -  1  )dX2  =  2en  dX2, 

which  implies  that  E  =  y/1  +  2cn  —  1.  For  small  strain  we  obtain  the  usual  relation, 

E  =  €n. 
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3.3.2  Global  Deformations  and  Elimination  of  Axial  Dynamics 

Consider  a  one  dimensional  structural  appendage  and  for  simplicity  consider  it  to  be  modeled 
as  a  beam.  Recall  that  each  infinitesimal  beam  element  is  identified  by  its  location  (i.e.,  the 
z-coordinate)  in  the  undeformed  configuration.  Thus,  the  element  configuration  variables 
rj(z,t)  and  £{z,t)  denote  the  spatial  location  and  rotation  at  time  t  of  the  infinitesimal 
element  located  at  coordinate  z  in  the  undeformed  configuration.  In  many  cases,  the  axial 
dynamics  (i.e.,  the  material  motion  along  the  beam  centerline)  evolves  on  a  substantially 
shorter  time  scale  than  does  the  transverse  or  shear  dynamics.  When  this  is  the  case  it 
is  often  desireable  to  reduce  model  complexity  by  eliminating  the  axial  dynamics  from  the 
equations  of  motion. 

One  approach  to  effect  this  reduction  is  to  enforce  additional  kinematic  constraints  on  the 
configuration  continuum.  The  approximation,  773  =  z  (cf.  Fig.  3.3)  is  a  commonly  employed 
approximation  and  valid  for  small  global  deformations.  Two  undesireable  consequences  have 
been  observed  as  resulting  from  this  approximation.  The  first  consequence  is  the  elimination 
of  foreshortening  effects  in  flexible  links  and  has  been  observed  to  be  of  importance  in 
precision  robot  placement  control  applications.  The  second  consequence,  which  may  be  of 
significance  in  large  angle  attitude  motions  of  flexible  space  structures,  results  from  the 
elimination  of  certain  inertial  coupling  effects.  It  is  important  to  emphasize  that  with  very 
large  structures  that  small  local  deformations  (i.e.,  strains)  may  result  in  significantly  large 
global  deformations  (i.e.,  displacements). 

A  variety  of  alternative  approximations  can  be  developed  which  eliminate  axial  dynamics 
while  retaining  correct  geometrical  relationships  during  structural  flexure.  We  explain  one 
such  approach  by  refering  to  Figure  3.3.  The  appendage  centerline  in  the  undeformed  state 
is  coincident  with  the  z  axis  and  has  length  L.  The  curve  CL  shown  is  the  centerline  of  the 
deformed  beam.  Notice  that  the  beam  centerline  at  any  time  t  is  parametrically  characterized 
by  the  map  CL  :  3?  — *■  SR3  given  by, 

V CL  ■=  V (M)> 


with  z  €  [0,L]  and  t  fixed.  At  any  point  z,  the  tangent  vector  to  the  curve  is 

drj 


If  we  insist  that  the  local  deformations  in  the  direction  of  vz  vanish  along  CL  then  we 
must  have 

ds2  -  ds20  =  0  (3.32) 

which  implies 

feMSMSr-* 


(3.33) 


Equation  (3.33)  should  be  thought  of  as  a  partial  differential  equation  in  773  which  is  to  be 
solved  subject  to  given  functions  of  z.  Notice  that  (3.33)  defines  a  unit  sphere  in 
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the  coordinates,  and  If  ^  is  assumed  to  be  positive,  then  we  may  take 


(3.34) 

and  the  solution  can  be  obtained  by  quadrature. 

3.4  Generic  Models  for  Slewing  and  Pointing  of  Precision  Optical  Structures 

In  this  section  we  develop  several  benchmark  generic  models  for  rapid  slewing  and  precision 
pointing  of  flexible  space  structures  which  are  motivated  by  problems  relating  to  control 
of  precision  optical  structures  subject  to  elastic  interactions.  Such  problems  arise  in  re¬ 
quirements  for  rapid  retargeting  coupled  with  precision  pointing  for  space-based  laser  (SBL) 
systems.  The  models  reflect  generic  qualitative  dynamical  properties  of  such  systems.  In  a 
subsequent  section  we  develop  a  simulation  model  with  physical  parameters  obtained  from 
the  benchmark  SBL  structural  model  developed  in  [Le87]. 

The  models  developed  in  this  section  focus  on  primary  sources  of  structural  interaction 
with  principal  body  slewing  maneuvers  affecting  system  LOS  pointing.  Modeling  assun  o- 
tions  used  to  characterize  generic  responses  are  based  on  the  initial  system  level  tradeoff- 
descril-’d  in  the  R  &  D  Associates  report  [Le87j.  This  study  indicates  that  the  principal 
source  of  structural  flexure  affecting  laser  LOS  is  withui  the  beem  expanded  optical  train — 
the  principal  structural  component  being  the  metering  truss  supporting  the  relative  position 
and  orientation  of  the  primary  and  secondary  mirrors.  Our  initial  or  first-level  model  assumes 
the  beam  expander  primary  mirror  and  support  is  rigidly  attached  to  the  spacecraft  body 
and  only  the  metering  truss  is  subject  to  flexure.  In  the  second  model  we  include  provisions 
for  articulation  of  the  SBL  beam  expander  with  respect  to  the  SBL  system  spacecraft  body 
using  a  gimbaled  joint. 


drj3 

dz 


L  Ah 


.  dr): 


./I  _  (^21)2  _ 

V  {dz  ’  Kdz} 
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3.4.1  Example:  Rigid  Body  With  1-Dimensional  Appendage 

We  consider  a  single  rigid  body  attached  to  a  flexible  appendage  as  illustrated  in  Figure  3.4. 
The  system  kinetic  energy  can  be  expressed  as 

T  =  \\Qu\\2 +  2(nu,  LR) +  \\R\\2 dm  (3.35) 

J  rb 

+  1  f  ||fiu  4-  u, ||2  +  2(f2u  +  ut,  LR)  +  ||i?||2<fm 

J  fb 

=  $u]£lbub  -  mbcTClbLR  + 

+  |  f  {11^6^  +  7t|| 2  +  2{CibTj  -f  1 jt,  LR)  4-  j|i?||2}/xd  dz 
J  c 

+  2/  +  £t)T I[ub  4-  £t]}p  dr 

Jo 

where  c  £  5f3  is  the  location  of  the  rigid  body  center  of  mass  in  the  body  frame,  mb  is  the 
mass  of  the  rigid  body,  Ib  is  the  inertia  tensor  of  the  rigid  body  in  the  body  frame,  and 
p  is  the  mass  density  of  the  beam,  rj(z,t )  is  the  position  vector  of  points  on  the  deformed 
centerline  of  the  beam  in  the  primary  body  coordinate  frame  and  £(z,t)  =  [ip,  9,  d>]T{z,  t)  is 
the  beam  angular  deformation.  We  have  assumed  small  deformation  of  the  beam  so  that  the 
angular  velocity  of  the  beam  section  at  z  is 

u(z,t)  =  wb  +  tt(z,t) 

up  to  first  order  in  the  angular  deformation. 

The  potentia1  energy  of  the  system  consists  only  of  the  potential  energy  associated  with 
deformation  of  the  beam.  Under  Timoshenko  beam  assumptions  [CKEFKPB68]  the  poten¬ 
tial  energy  function  is 

v((,v)  =  k£{aj(*.r+Ew#+Ei*+.r 

4-  k\GA(t]iiZ  —  9 )2  -|-  K2GA(ri2>z  4"  <f>)2  4"  A(Tj3tZ  —  l)2}  dz  (3.36) 
=  k  f‘  {(£**.)  +  Ofr  -  Pt)TS(v,  -  Pi)}  dz 

where  the  stiffness  matrices  are  defined  as 

K  =  diag(G7,  EI2,  EI3),  S  =  diag(/tiGA,  K2GA,pEA), 

and 

‘  0  1  0  ‘ 

P=  00-1  . 

.00  0 

The  system  angular  momentum  vector  is 

a  =  Ibui  4-  f  [At]  x  (ub  x  T}  4  7)  4-  Iv]p  dz 
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so  that  equations  (3.19)-(3.21)  reduce  to 
rr(£)(/&ci>fc  +  u>b  x  Ib^b) 

+  I  [rT(0^  x  ivtt  +  u>b  X  (&b  X  Tf)  +  u>b  X  T]  +  2wb  X  Tj) 

Jo 

+  Atj  x  LR  +  I(uj  +  u>b  x  u>)\p  dz  -f  me  x  LR  =  Q(b,  (3.37) 

mR  +  mLT (u)  x(wxc)-|-«xc)  +  /  pA(LT D2tj  -f  R)dz  =  Qr,  (3.38) 

Jo 

pA(rjtt  +  u>b  x  {ujb  x  tj)  +  ujb  x  jj  -f  2cob  x  77  +  LR)  =  Qv  -  5(7/zz  -  P£z),  (3.39) 

pl[itt  +  Ub  +  Ub  X  ( u>b  +  6)]  =  Q(-  Kizz  +  PTS(tjz  -  P£).  (3.40) 

System  Dissipation  A  simple  model  of  generalized  dissipation  in  the  appendage  can  be 
obtained  by  introducing  a  Rayleigh  dissipation  function.  We  formulate  such  a  function  based 
on  the  assumption  that  dissipation  forces  are  proportional  to  beam  deformation  rates,  i.e. 
generalized  coordinate  velocities  ( Tjt  and  (t)  and  strain  rates  ((771  )tz  —  0t,(rj2)tz  —  <f>t,(tz)- 

R(Vt,tt)  =  k  [  {vfzilt  +  +  (v tzf^aVtz  +  (6*)rH4(6*)}dz  (3.41) 

Jo 

where  E;  =  diag(^,  fa,  £3).  From  R(rjt,^t)  we  obtain  the  generalized  dissipative  forces 


—  =  ElVt  +  E3(Vtzz-P^tz), 

OVt 

SR  „  _  .  ne  \  — 

77"  =  ci2  +  z.3\Vtz  —  Pit)  +  ^iStzz  ■ 


(3.42) 

(3.43) 


3.4.2  Example:  Articulated  Bodies  With  Flexible  Appendage 

We  now  consider  a  modification  of  the  previous  Example  which  includes  a  second  body 
attached  to  the  primary  body  with  a  three  axis  gimbal  as  illustrated  in  Figure/2)  In  addition, 
the  second  body  carries  with  it  a  reaction  wheel  package.  The  kinetic  energy  function  is 

T  =  ^J/pu;p  -  mpCpUpLR  +  £mp||P||2 

' - v - - 

primary  body 

-I-  iwj  I,u>,  -  m,cT,OtLt{LR  -|-  flpCffm}  -f  \m,\\LR  +  OpCpm||2 

' - v - ' 

second  body  (wheels  locked) 

+  I3UJ3  +  U)J  I3U)3 

' - v - ' 

reaction  wheels 


+  !  J‘{\\0p<l  +  Vtf  +  2(!5PV  +  1.,  LR)  +  iz 

+  S  j  {("»  +  (i)Tj(up  +  ( 1 )}  Pdz 


(3.44) 
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Figure  3.4:  Rigid  Body  with  Flexible  Appendage 


besm  exponder 


Figure  3.5:  Hinged  Bodies  with  Flexible  Appendage 
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Note  that  cp  and  c,  are  the  locations  of  the  center  of  gravity  of  *he  primary  body  and 
second  body  (ineluding  wheels)  in  their  respective  body  coordinate  frames.  The  vector  cgm 
denotes  the  location  of  the  three  axis  gimbal  in  the  primary  body  frame.  The  matrix  I\ 
is  the  primary  body  inertia  tensor  in  the  primary  body  frame,  I2  denotes  the  second  body 
inertia  tensor  with  reaction  wheels,  and  /3  is  the  diagonal  matrix  of  wheel  inertias. 
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4  Feedback  Linearization  of  Lagrangian  Systems 

4.1  Output  normal  form  equations  for  Lagrangian  Systems 
Consider  the  case  of  a  finite  dimensional  Lagrangian  system  with  Lagrangian, 

L{q,q):M  xTqM  »  (4.1) 

and  assume  that  the  configuration  manifold,  M.  is  locally  diffeomorphic  to  9?n.  Thus  the 
n-vector  q  has  elements  which  are  the  generalized  coordinates. 

For  control  design  we  assume  there  are  a  finite  number  of  external  generalized  forces 
which  can  be  exerted  on  the  system  which  we  collect  in  the  m- vector  /  6  3f?m.  The  external 
forces  are  introduced  by  the  definition  of  a  differential  virtual  work  expression  in  terms  of 
the  generalized  coordinates;  8W  =  £,■  Qfdqi.  We  assume  that  the  external  forces  in  physical 
coordinates  are  such  that  m  <  n  and  for  each  coordinate,  Qi  =  =  Gif  and  the  virtual 

work  takes  the  form, 

SW  =  J2  fGfdqi  =  fTGTdq,  (4.2) 

t=l 

where  GT  =  [Gf , . . . ,  G£]. 

Assumption  1:  The  LagTangian  system  is  regular  with  respect  to  external  generalized 
forces;  i.e.,  the  operator  G  is  rank  m  in  a  neighborhood  of  every  configuration  q  £  M. 

Next  we  identify  a  set  of  system  primary  outputs,  y  6  9?n,  given  by, 

V  =  M$)»  (4-3) 

where  h  :  M  —*  is  a  G°°  function.  We  also  assume  that  there  exists  a  C°°  function, 
t(q)  i-*  u,  such  that  U(q)  =  h(q)®t(q)  is  a  local  diffeomorphism;  i.e.,  there  exists  an  inverse 
map; 

q  =  g{y,u)  (4.4) 

g  :  x  Mu  -*  M. 

Proposition  1:  Under  the  above  assumptions  the  Lagrangian  system  can  be  decomposed 
and  written  in  explicit  output  coordinates  form 7; 

(4.5) 

(4.6) 
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where 


L(y,y,u,u) 

9(u,y) 


9)l9=g(ti,v),9=g(u,v)  » 

Dy(g)y  +  Du(g)u. 


(4.7) 

(4.8) 


Proof  of  Prop  1:  To  see  this  note  that  (4.8)  follows  from, 

d  [<9<7]  .  \dg]  . 

Also,  the  virtual  work  expression  is  rewritten  in  the  new  coordinates  as, 

SW  =  fTGTW\iv  +  fTGT\^\iu. 
ay  ou 

Then  direct  application  of  Lagrange’s  equations  in  the  ( y,u )  coordinates  gives  the  desired 
result. 


□ 


4.2  Structural  Assumptions  on  Lagrangian  Systems  and  PLF  Compensation 

From  the  above  section  we  see  that  under  fairly  general  assumptions  a  Lagrangian  system 
with  n  degrees  of  freedom,  m  external  forces  and  m  outputs  defined  as  in  (4.3),  that  the 
resulting  equations  of  motion  can  be  obtained  by  application  of  Lagrange’s  equations  in  the 
explicit  output  form  (4.5)-(4.6)  which  we  write  as, 

Jy{y,u)y  +  N(y,u)u  +  <f>v(y,  y,u,u)  =  Gyf  (4.9) 

Ju(y,u)u  +  NT(y,u)y  +  <f>u(y,y,u,ii)  =  Guf.  (4.10) 

Assumption  2:  The  system  inertias;  i.e,  Ju,  Jy, 

nite  operators  for  q  G  M. 

Then  we  have  Jy  —  NJ~*NT  is  also  positive  definite  on  M.. 

Definition  1:  We  say  that  the  above  m-input  by  m-output  system  has  local  relative  degree 
2  for  the  configuration  (yo,u0)  if  and  only  if  the  m  x  m  matrix,  [Gv  -  N J~lGu){y,u)  is 
invertible  everywhere  in  a  neighborhood  of  (y0,«o)- 

A  Partial  Linearizing  Feedback  (PLF)  Compensation  is  a  nonlinear  feedback  transforma¬ 
tion  of  the  form, 

/  =  £v(y,  y, «.  u)  +  £a(y,  y , u)a 

with  Ky,Ka  :  M.  x  TqM  — ►  J?m,  which  renders  the  compensated  response  in  y  linear  and 
time  invariant  in  a,  and  decouples  the  remaining  (zero)  dynamics  (u,  u). 


Jy  N 
NT  Ju 


are  symmetric,  positive  defi- 
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Proposition  2:  Assume  the  Lagrangian  system  has  local  relative  degree  2  for  every  con¬ 
figuration  q  €  A  C  M  where  A  is  an  open  dense  subset  of  M .  Then  a  Partial  Linearizing 


Feedback  (PLF)  compensation  is  a  feedback  of  the  form, 

/  =  -[(?„  -  iVJ-'GJ-1  {*„  -  +  [J,  -  NJ;'NT)a}  (4.11) 

transforms  the  closed  loop  system  to  input/output  linearized  and  decoupled  form; 

y  =  a,  (4.12) 

Ju{y,u)u  +  N(y,u)Tu  +  <f>u{y,y,u,u)  =  Gvua  (4.13) 

where 

4>l  :=  <t>u  +  Gu[Gy  —  N J~* Gu]'1  (<j>y  —  NJ~l(j)u),  (4.14) 

GZ  :=  Gu[Gy  -  NJ~xGu\~l[NJ~lNT  -  Jy\,  (4.15) 

where  a  €  represents  the  control  inputs  in  “acceleration  coordinates”. 


The  system  primary  outputs  y  given  by  (4.3)  induce  a  decomposition  of  the  system 
configuration  space  as  M  =  My  x  Mu  where  dim  My  =  m  and  Mu  =  {q  (z  M  :  h(q)  =  0}. 
The  PLF  compensation  (4.11)  decomposes  the  trajectories  of  (4.9)-(4.10)  into  controlled 
responses  on  My  x  TvMy  and  output  constrained  or  zero  dynamics  on  Mu  x  TUMU.  Then 
following  Byrnes  and  Isidori  [BI85]  we  say  that  the  system  (4.9)-(4.10)  is  locally  minimum 
phase  if  the  output  constrained  dynamics  are  asymptotically  stable  to  the  origin  ( u,u )  = 
(0,0). 


Proposition  3:  The  system  (4.9)-(4.10)  is  locally  minimum  phase  if  and  only  if  the  system, 


Jj(a)ti  +  ^(»,i)  =  0 


(4.16) 


is  asymptotically  stable  to  the  origin,  ( u ,  it)  =  (0, 0),  where 


^u(y,y,«,«)lv=0,i=o -  -NJU'GU]  x(<j>y 

J«(y>u)\v=o,v=o  • 


y=0,y=0  ’ 


(4.17) 

(4.18) 


Proof  of  Prop’s  2  &  3:  To  see  this  result  we  first  solve  for  the  accelerations  in  the  pair 
of  equations  (4.9)-(4.10)  to  obtain, 


(4.19) 


40 


Techno-Sciences,  Inc. 


TSI-90-10-10-WB 


where  we  have  assumed  that  the  system  inertia  is  symmetric,  positive  definite  everywhere 
in  ( y,u )  €  My  x  M.u.  The  desired  decoupling  and  linearizing  control  is  then  given  as, 


/  =  r;'(-s,  +  a) 

Now  recall  that  under  the  assumption  on  the  system  inertias  we  have, 


'  Jv 

N 

-1 

A"1 

nt 

j u  . 

-J;lNT  A'1  j-1  +  J~1NtA~1NJ~1 

where  A  =  Jy  —  NJU  1NT.  Thus  we  find 

S,  =  A"Vb  -  A-^J-Vu, 

rv  =  A-1  -  A^NJ^Gu. 

Furthermore,  since  we  have 

r;%  =  [a-‘g,  -  a-'jyj^gj-^a-1*,  - 
=  [Gy  -  NJ;'Gy]-'(<h  -  NJ?*y), 


(4.20) 

(4.21) 

(4.22) 

(4.23) 

(4.24) 


then  (4.11)  follows  from  (4.20). 

To  show  (4.12)-(4.15)  substitute  (4.11)  into  (4.10).  Likewise,  (4.16)-(4.18)  follow  from 
(4.12)-(4.15)  by  setting  a  =  0  which  from  (4.12)  implies  that  y  =  y  =  0. 


□ 


Remarks:  The  objective  of  PLF  compensation  (4.11)  is  to  obtain  input/output  lineariza¬ 
tion  of  a  y  by  decoupling  of  the  u  dynamics  from  the  outputs,  y.  We  assert  that  for  the 
Lagrangian  systems  considered,  under  the  assumptions  1  and  2,  and  assuming  that  for  the 
given  /  and  y  variables  that  the  system  is  locally  minimum  phase,  then  the  response  of  the 
PLF  compensated  system  consists  of  a-controlled  motions  on  Mv  and  stable  transient  mo¬ 
tion  converging  asymptotically  to  the  origin  on  Mu.  For  regulation  of  the  primary  outputs 
to  zero  the  dynamic  response  is  decomposed  to  asymptotic  motions  toward  the  manifold  My 
described  by  (4.12)  and  motions  on  My  given  by  (4.16). 

Design  of  the  a  control  is  of  course  simplified  by  the  linear  dynamics  given  by  (4,12). 
As  most  flexible  structure  systems  are  Lagrangian  systems  with  local  relative  degree  two, 
the  design  in  acceleration  coordinates  is  particularly  attractive  feature  of  the  method.  The 
PLF  compensation  (4.11)  implements  an  online  Inverse  Force  Model  (IFM)  which  effects  the 
transformation  to  acceleration  coordinates.  In  the  sequel  we  show  that  the  stable  implemen¬ 
tation  of  the  IFM  may  be  sensitive  to  parasitics  in  the  decoupled  dynamics  (4.13). 
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5  Reduced  Order  Models  and  PFL  Control  of  Flexible  Struc¬ 
tures 

Models  obtained  by  finite  dimensional  approximation  of  the  Euler- Lagrange  equations  or  of 
the  system  Hamiltonian  as  outlined  above  will  typically  obtain  a  large  order  finite  element 
model  in  the  form, 

M{q)q+ B{q,q)q+ K(q,q)q  =  Gf,  (5.1) 

where  we  assume  the  system  inertia,  M(q)  :  M.  — *  3inXn  is  symmetric,  positive  definite  on 
M,  reduced  configuration  space  with  dimension  n  and  stiffness  K(q,q )  :  M  x  TqM  — ►  §?nXn 
is  symmetric,  positive  semidefinite  stiffness,  and  B(q,  q)  :  M  x  TqM  — ►  9?nxn  models  system 
dissipation  on  M.  x  TqAi. 

Our  interest  in  model  order  reduction  arises  from  consideration  for  robust  implementation 
of  PLF  compensation  for  flexible  structures.  In  general,  the  dynamics  are  nonlinear  and 
we  base  reduction  on  time  scale  decomposition  with  respect  to  a  minimum  energy  state 
of  the  mechanical  structure  model.  For  space  structure  applications  there  are  typically 
cyclic  coordinates  involving  attitude  rotations  which  do  not  appear  explicitly  in  the  system 
Lagrangian  [Gol82].  These  coordinates  can  be  translated  to  the  origin.  Thus  we  focus 
attention  on  the  linear  perturbation  model  about  the  origin  in  M  x  TqM  to  obtain  time 
scale  decomposition. 

Fact:  Under  the  above  assumptions  there  exists  a  nonsingular  n  x  n  matrix,  P,  such  that 
PT MqP  =  PT KqP  =  Kd  =  diag{o/j, . . .  ,u>*} 

with  u>i  <  •  ■  •  <  ujn,  where  M0  =  M{ 0),  K0  =  /f(0, 0). 

Assumption  3:  Given  P  as  above,  let  B0  =  5(0, 0)  and  assume  that8 

PT B0 P  =  Bd  =  diag{2(1u;1, . . . ,  2<nu’n}- 


The  transformation  q  —  Px  obtains  the  linear  perturbation  model  in  decoupled  “modal” 
coordinates  [Bal78], 

x  +  Bdx  +  Kdx  =  Gf,  (5.2) 

where  G  =  PTG,  whereas,  the  original  nonlinear  equations  (5.1),  written  in  modal  coordi¬ 
nates: 

M{x)x  +  B(x,x)x  4-  K(x,x)x  =  Gf,  (5.3) 

8The  assumption  of  “modal  damping”  is — although  typical  in  modeling  flexible  space  structure 
dynamics — an  assertion  which  may  require  careful  validation  in  practice.  Among  other  limitations  this 
assumption  implies  that  gyroscopic  effects  are  negligible. 
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where 

M(x)  =  PT  [M(q)\^Pl]  P,  (5.4) 

K(x,i)  =  Pr[jr(t,«)|„*,irf4]f',  (5.5) 

B(x,i)  =  PT  [S(9,9)|q=Pllj=Pi]  P,  (5.6) 

may  be  highly  coupled  in  a  neighborhood  of  the  origin. 

The  standard  approach  to  model  order  reduction  for  linear  models  of  flexible  structures 
is  to  retain  the  first  p  modes  and  neglect  the  n  —  p  residual  states.  Such  an  approximation 
assumes  that  the  system  modal  coordinates  are  decomposed  as  x  =  [(x°)r,  (x1)7]1^  where  x1 
is  a  fast  transient  (or  parasitic)  in  the  time  scale  of  x°.  Thus  we  seek  to  approximate  the 
system  potential  energy  by  the  decomposition  of  the  stiffness  in  (local)  modal  coordinates 
model  (5.2); 

Kd  = 

To  transform  the  model  (5.2)  into  the  standard  form  for  singular  perturbation  analysis 
we  must  identify  an  appropriate  scaling  of  the  stiffness.  One  natural  choice  of  scaling  is  to 
let,  Kd  —  ^Kd°,  and  the  residual  coordinates,  x1  =  e2z,  such  that  z  ~  0(1)  in  e2.  For 
example,  we  may  take, 

€(p)  =  — —  (5.7) 

wp+i 

for  which  we  see  that  e  —*  0  as  p  is  increased.  Then  by  assumption  3  we  have,  Bd  =  7 5™, 
B°  0 

where  Bd  =  d  ;  a  fact  which  will  become  useful  in  the  sequel. 

0  “d 

Then  (5.2)  is  decomposed  into  p-dimensional  retained  coordinates,  x°  and  (n  —  p)- 


dimensional  residual  coordinates,  x1, 

«  =  [ft, -P.l  (*”.),  (5-8) 

where  P0  :  »p  ->  3T,  Pl  :  »n“p  -►  ft",  such  that 

x°  +  B°dx*  +  K°x°  =  Go/,  (5.9) 

x1  +  B^x1  +  K\xl  =  Gif,  (5.10) 

which  after  scaling  becomes  singularly  perturbed  [ABB89]  as, 

x°  +  Bdi°  +  Kdx°  =  Go/,  (5.11) 

e2[z  +  B\z\  +  K™  z  =  Gif.  (5.12) 

After  decomposition  and  scaling  the  FEM  model  (5.3)  is  also  singularly  perturbed, 

Moo*°  +  e2Moi*  +  Boo*0  +  e2Boii  +  Aoo*°  +  e2Koiz  =  Go/  (5.13) 

M«x0  +  e2Mnz  +  B10*°  4-  +  K^x°  +  K™z  =  Gif.  (5.14) 
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The  reduced  order  model  is  obtained  as  the  quasi-steady  state  or  slow  time  system  by 
letting  €  =  0; 


where 


+  B°mx°  +  K^x- 

=  G0f, 

(5.15) 

[AC|V  +  B°„x°  +  «|V  +  K\°z 

=  Gif, 

(5.16) 

AC(*“)  = 

(5.17) 

=  P?B(q,q)P„ 

\q=P0x°  ’ 

(5.18) 

kS,(x°,i°)  =  P*K(q,q)P„  |a=iW, 

A=Po±°  ’ 

(5.19) 

Qi 

O 

II 

P 

(5.20) 

where  (5.15)  is  the  reduced  model  [Kok87j.  By  assumption,  K™  is  invertible  and  we  can  solve 
explicitly  for  the  slow  time,  or  quasi-steady  state  component  of  z; 

z  =  \K™]~x{Gf  -  Bx°  -  Kx0},  (5.21) 

where 

B  =  Bto-M&Mn'Boo,  (5.22) 

K  =  K&  -  M^MooKoo,  (5.23) 

G  =  Gi-  Ml  Moo1  Go-  (5.24) 


The  boundary  layer  correction;  z  =  z  +  77  is  obtained  in  the  fast  (stretched)  time  scale, 
r  =  t/e  where,  ^  1  ^  and  letting  e  =  0  after  substitution  in  (5. 13)— (5. 14).  Introducing 

the  notation,  x'  :=  we  can  write  the  fast  time  scale  system  as, 

Muv"  +  Bl0lV'  +  K™n  =  Gif.  (5.25) 


A  standard  result  in  singular  perturbation  theory  is  that  the  reduced  model  (5.15)  is  0(1) 
approximation  for  e  C  1  if  the  fast  time  scale  system  (5.25)  is  asymptotically  stable  [KKO86]. 
Then  the  time  scale  decomposition  suggests  that  trajectories  of  (5.1)  consist  of  fast,  transients 
of  (5.25)  which  converge  to  the  manifold  So  =  {(*,  x)  G  Mo  x  TxMo  :  z  =  z,  z  =  0},  where 
dim  A4o  =  V  and  trajectories  of  (5.15)  on  So- 


Remarks:  For  flexible  space  structures  we  typically  assume  that  the  stability  of  the  fast 
time  scale  model  is  due  to  inherent  damping  of  the  structure  rather  than  the  introduction 
of  rapid  control.  One  issue  of  practical  concern  is  the  controllability  of  (5.25),  which  for 
the  case  of  flexible  structures  depends  on  the  choice  of  p  and  slowly  maturing  technology  for 
wide  band  vibration  control.  When  the  fast  time  scale  system  is  only  weakly  controllable  (i.e., 
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controllable  only  through  the  slow  time  system)  the  contribution  of  these  modes  is  clearly 
“parasitic”  and  neglecting  these  dynamics  will  contribute  to  the  robustness  of  the  control 
design  [Kok87].  Model  order  reduction  in  flexible  space  structure  control  design  is  typically 
a  multistep  process.  We  recommend  that  weak  controllability  issues  related  to  actuation 
technology  be  considered  as  part  of  any  initial  reductions.  Model  order  reduction  for  control 
design  for  flexible  structures  is  thus  effected  by  the  technology  available  for  actuation  and 
sensing.  Structural  input/output  properties  which  also  effect  choice  of  model  order  reduction 
methods  include  weak  vs.  strong  controllability  and/or  observability  of  the  reduced  models. 
These  and  many  other  issues  can  effect  the  choice  of  time  scaling  and  model  reduction. 

5.1  PLF  Compensation  for  Reduced  Models 

Reduced  order  models  of  varying  dimension  can  be  obtained  by  the  time  scale  analysis  of 
the  previous  section  where  the  approximate  system  dynamics  evolve  on  the  reduced  man¬ 
ifold  £o  C  Mo  x  TqM0-  In  the  next  few  paragraphs  we  summarize  the  transformation 
of  the  reduced  model  to  output  normal  form  and  the  construction  of  PLF  compensation 
based  on  the  reduced  model.  Our  objective  is  to  demonstrate  how  considerations  for  robust 
implementation  of  PLF  compensation  for  flexible  structures  will  drive  the  choice  of  model 
reduction. 

We  focus  attention  on  an  m- vector  of  system  outputs  (4.3)  for  the  Lagrangian  system  to 
the  form  (4.5)-(4.6).  The  output  map  is  given  in  reduced  model  coordinates  by9, 

y  =  h\x)  :=  M?)U0X  • 

Then  the  inverse  map  for  the  reduced  order  model  is  from  (4.4), 

x  =  [P0r  p0]_1-P0r9(y»«)  =:  9°{y,u)- 

Under  smoothness  assumptions,  we  can  write, 

9°{y,  u)  =  y°(y,  u)y  +  g°u(y,  u)u,  (5.26) 

and  following  the  procedure  of  §  5  we  identify  the  reduced  model  (5.15)  in  the  form  (4.9)- 
(4.10)  where, 

(5.27) 

(5.28) 

(5.29) 

9In  this  section  we  replace  the  slow  time  state  x°  €  J?p  with  x  for  notational  convenience. 


Jy  -  ( 9l)TMgl,  Ju  =  [gl}TMg°u,  N  =  [g°y}TMgl 
Gv  =  [g°v}TG,  Gu  =  [g°u]TG, 

'4>V  \  =  \  Byy  Byu  1  A/  \  [  Kyy  Kyu  1  (V  \ 

k4>u)  BUy  Buu  j  \U)  [  KUy  Kuu  J 
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and 


Byy 

Byu 

Buy 

BUu 

Kyy 

Kyu 

Kuy 

Kuu 

(y,y,u,u) 


(l t,V,u,u) 


w- 1 


01r 


LbS 
i  Wf 


K[g°y,g°u  ]. 


(5.30) 

(5.31) 


The  PLF  control  on  Aio  x  TxMo  is  given  by  (4.11)  with  (5.29).  The  analysis  of  section  2.2 
carries  over  directly  to  the  decomposition  of  the  reduced  manifold,  50- 

Let  Mq  C  M  be  the  p-dimensional  reduced  configuration  space  obtained  by  retaining 
the  p-dimensional  slow  time  modes  according  to  the  scaling  e(p)  as  in  (5.7).  Then  given  a  se¬ 
quence  of  integers,  pk,  such  that,  n  >  px  >  •  •  •  >  p*.  >  m  the  reduced  system  manifold  satisfy, 
M.  C  Mok  C  and  .Ado*  =  Adyk  ©  where  dim  Mvyk  =  m  and  dim  M„k  =  pk  —  m.  In 
this  sense  we  understand  that  when  p  =  mwe  have  the  case  usually  considered  in  robotics 
[SV87]  where  the  nominal  (reduced)  model  is  exactly  feedback  linearizable  by  nonlinear 
transformation  of  the  system  state  together  with  feedback.  For  these  cases  a  variety  of  con¬ 
trol  design  techniques  have  been  devised  to  assure  robust  stabilization  subject  to  bounded 
model  uncertainty  of  the  nominal  plant  model  described  by  a  regular  perturbation  [Spo87]. 
Unlike  the  robotic  problems,  flexible  structure  control  problems  are  characterized  by  high 
dimensional  reduced  order  models  where  m  <  n  typically.  The  IFM  implemented  by  PLF 
compensation  obtains  linearization  through  transformation  of  the  design  model  to  accelera¬ 
tion  coordinates.  A  critical  feature  for  practiced  implementation  of  PLF  control  for  flexible 
structures  is  the  tradeoff  between  the  fidelity  of  the  online  IFM  vs.  the  model  sensitivity. 
We  will  highlight  these  issues  by  reference  to  the  case  when  the  zero  dynamics  of  the  high 
order  design  model  are  singularly  perturbed. 


5.2  Model  Reduction  for  Robust  Implementation  of  PLF  Control 

One  way  that  nonminimum  phase  responses  can  arise  in  flexible  structure  systems  is  due 
to  noncollocation  of  localized  point  forces  and  localized  position  outputs.  It  is  apparent 
that  model  order  reduction  can  in  some  cases  reduce  a  nonminimum  phase  system  to  one  of 
reduced  order  which  is  minimum  phase.  This  is  significant  for  robust  control  design.  Often 
instabilities  of  the  zero  dynamics  are  associated  with  fast  time  scale  effects  which  may  be 
safely  neglected.  Indeed,  as  shown  by  Sastry  et  al  [SHK89]  regular  perturbations  of  the 
system  model  for  both  linear  and  smooth  nonlinear  systems,  affine  in  control  can  induce 
singular  perturbations  of  their  zero  dynamics.  Thus  to  determine  whether  the  fast  time 
transients  of  the  zero  dynamics  for  a  given  reduced  model  are  stable  or  not  may  require 
highly  accurate  system  models.  It  is  clear  then  that  good  engineering  practice  and  robust 
design  principals  will  suggest  that  fast  time  transients  of  the  zero  dynamics  for  a  given 
reduced  order  model  should  be  reduced  from  the  system  model  before  PLF  implementation. 
Thus  we  recommend  that  the  system  model  be  further  reduced  when  the  zero  dynamics  exhibit 
time  scale  separation. 
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The  zero  dynamics  (4.16)  can  be  time  scaled  if  there  exists  a  transformation  on  Mu  of 
the  form  (^1 ,  ^2 )  =  t(u )  such  that  the  model  is  decomposed, 

+  JlJ,  +  =  0,  (5.32) 

W./«i  +  JlJ*  +  <,  =  0,  (5.33) 

and  can  be  scaled,  £2  =  ez  so  that  z  is  the  fast  time  state  and  the  model  is  singularly 

perturbed: 

+  eJu, 12*  +  ft,  +  efti  =  0,  (5.34) 

[^,12]  6  +  e^u,22^  +  ft2  4-  =  0.  (5.35) 

This  decomposes  the  zero  dynamics  manifold  as  Mu  =  M^  ©  M^  where  dim  M{x  =  p  —  m, 
dim  =  n  —  p. 

The  implication  for  reduced  order  PLF  control  is  now  considered.  Reduced  order  PLF 
compensation  can  now  be  designed  for  the  reduced  model  obtained  from  the  decomposition 
of  the  configuration  manifold,  M  =  My®M^x  ©  M(3  and  projection  onto  the  reduced  space, 
M  =  My  ©A4{,.  The  principal  question  which  arises  is  under  what  conditions  can  the  time 
scaling  of  the  zero  dynamics  lead  to  time  scale  decomposition  of  the  Lagranian  system  in 
which  the  fast  time  scale  system  can  be  neglected. 

Following  the  ideas  in  [B185]  we  assert  that  there  exists  a  coordinate  transformation  on 
Mu  such  that  the  explicit  output  coordinates  expression  (4.9)-(4.10)  is  in  output  normal 
form  and  under  the  previous  assertion  the  dynamics  on  Mu  can  be  decomposed  so  that  we 
can  write  the  system  dynamics  in  the  form, 

■Jy(y.0y  +  ^v(y»y»£»0  =  °y /»  (5.36) 

Ju,n(y,()(\  +  Ju,n(y,0^  +  </>u,i{y,y,L0  =  (5.37) 

[«/u.w(y,  0]r6  +  A32(y,06  +  =  0.  (5.38) 

We  assert  that  the  above  normal  form  equations  are  also  singularly  perturbed  and  the  re¬ 
duced  model,  obtained  by  projection  onto  the  configuration  space  A4q  =  My  x  M is  a 
valid  slow  time  scale  approximation  ets  e  — »•  0  if  the  boundary  layer  system  is  asymptotically 
stable.  Note  that  the  boundary  layer  system  for  the  normal  form  equations  will  have  dif¬ 
ferent  stability  properties  in  general  Rom  the  boundary  layer  system  for  the  perturbed,  zero 
dynamics  (5.34)-(5.35).  Indeed,  even  if  the  zero  dynamics  have  fast  time  responses  which 
are  unstable  then  we  assert  that  PLF  control  is  implementable  for  the  reduced  order  model 
obtained  by  projection  of  the  system  dynamics  onto  My  x  M^x  if  the  fast  time  scale  part  of 
(5.36)-(5.38)  is  asymptotically  stable. 

Remark:  The  model  order  reduction  obtained  by  time  scale  decomposition  of  the  zero 

dynamics  when  applied  to  the  system  normal  form  equations  obtains  a  p  DOF  reduced 
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model  which  matches  the  p  slow  time  modes  of  the  zero  dynamics.  Thus  the  system  relative 
degree  assumptions  will  be  unchanged  under  this  form  of  model  reduction.  The  reduced 
PLF  control  will  implement  a  slow  time  scale  approximation  of  the  IFM. 
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6  Experiments  in  Slewing  and  Pointing  of  Multibody  Systems 

Our  studies  during  the  course  of  this  project  have  investigated  some  of  the  potential  benefits 
that  can  accrue  from  implementation  of  dynamic  inversion  by  PLF  for  control  of  flexible 
structures.  On  the  surface,  achievable  performance  of  PLF  compensation  depend-’  on  the 
accuracy  of  models  used  to  predict  the  nonlinear  couplings.  However,  during  the  course 
of  this  project  we  have  focused  attention  on  the  pragmatic  issues  of  control  design  with 
model  uncertainty.  In  this  section  we  describe  the  assumptions  used  to  obtain  reduced 
models  for  multiaxis  slewing  and  pointing  control  of  an  SBL  system.  We  describe  several 
simulation  experiments  which  can  provide  a  basis  for  design  of  laboratory  experiments  to 
validate  the  predicted  performance  of  PLF  compensation  for  slewing  and  pointing.  We  have 
included  available  data  on  the  ASTREX  facility  and  the  first  test  article  in  our  modeling, 
control  design,  and  simulation  studies  in  an  effort  to  determine  initial  feasibility  of  conducting 
experiments  to  validate  nonlinear  control  laws  on  this  facility.  In  the  following  paragraphs  we 
describe  a  series  of  experiments  which  could  be  run  on  ASTREX  and  used  to  validate  SBL 
slewing  and  pointing  control  with  PLF  methods  used  to  compensate  for  dominant  system 
nonlinear  couplings. 

6.1  Multibody  System  Model  for  Multiaxis  Slewing  and  Precision  Alignment 

The  modeling  method  of  Lagrange’s  equations,  as  described  in  the  report  [BBKA88],  along 
with  spatial  discretization  via  collocation  by  splines  to  develop  a  finite  dimensional  model 
suitable  for  simulation  of  large  angle,  multiaxis  motion  of  a  generic,  two-body  model  of  a 
SBL  system  beam  expander  is  described  in  this  section.  Beginning  with  the  generic  flexible 
spacecraft  model  [BBKA88,  sec  4]  the  following  model  assumptions  are  made: 

1.  axial  appendage  deformations  are  negligible10,  773  «  z, 

2.  translation  velocity  of  the  system  primary  body  is  negligible,  R  &  0, 

3.  torsional  appendage  deformations  are  negligible,  rf>(z)  ~  0. 

Thus,  following  the  DPS  modeling  approach  discussed  in  [BBKA88,  sec  3],  the  configuration 
space  for  the  distributed  parameter  model  is  50(3)  x  Hg,  where  Hg  is  the  set  of  continuously 
differentiable  functions,  x{z)  —  [*7ii *7zi 0, <f>]T  £  Hg  defined  on  the  interval  z  6  [0,f]',  and 
which  satisfy  the  geometric  boundary  conditions  : 

Q  '■  ’/ll2)  =  0,  tj2 (z)  =  0,  0(z)  =  0,  <f>(z )  =  0,  at  z  =  0. 

With  the  notations!  conventions  in  [BBKA88],  the  model  will  be  described  by  the  appendage 
lateral  deformation,  f](t,z)  =  along  the  body-fixed  *  and  y  axes,  respectively,  and 

10See  §  3.3  for  a  discussion  of  this  assumption. 
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appendage  angular  deformation,  £( z,t )  =  [#(z),  <f>(z)\T .  The  attitude  angular  rates  referenced 
to  the  body  fixed  frame  are  denoted, 

The  potential  functions  for  the  variational  analysis  leading  to  the  dynamic  equations  of 
motion  reduce  in  this  case  to: 


T  =  IfyUlb  (6-1) 

+  k  [  fl^z)  +  Vt(z)\\2pA  dz 
J  0 

+  5  Jo  K  +  Ptt{z)]TI[u b  +  Pit{z))p  dz , 


where 


v(i»,o  =  if‘{Ei,el(z)  +  Ewl(z) 

+«iCM[(t7i(.z))i  -  6(z) ]2  +  k2GA[(tj2(z))z  -  4>{z)]2}  dz , 

R{*hi)  =  £  Jo  {vI^TJt  +  (fz2tt  +  ( Vtz)TZ3T]tz  +  (ttz)TZ4(Ztz)}dz, 


ro  0  n 


p  = 


0  1  0 
1  0  0 


(6.2) 

(6.3) 


T  is  the  kinetic  energy,  V  the  potential  energy,  and  R  the  Rayleigh  dissipation  function.  The 
notation  used  is  standard  in  continuum  mechanics  of  beams  and  is  summarized  in  Table  6.1. 
The  area  moment  tensor  is  assumed  to  have  the  form11, 


'(*) 


'/**  0  0  ' 

0  Iyy  0 

.  0  0  Izz  . 


For  simulation  studies  contained  in  this  report  we  utilize  a  simple  damping  model  given 
by  the  assumption  that  the  matrix  coefficients  in  the  dissipation  function  are  of  the  form, 

-  _  Cu  0  ' 

“*  “  0  <„•  • 

Note  that  Si, £2  model  external  dissipation  while  £3,54  model  internal  dissipation  [BK89). 
The  control  forces  acting  on  the  system  primary  body  are  modeled  as  external  torques  and 
enter  the  variational  setup  through  a  virtual  work  expression  of  the  form, 

SW  =  rf6ih.  (6.4) 

11  We  use  the  so-called  NASA  standard  or  321  convention  [Gol82]. 
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Notation 

Explanation 

P 

mass  density 

A 

cross  section  area 

E 

elasticity 

kG 

effective  shear  modulus 

I 

area  moment  of  inertia 

•  dx 

x  —  dt 

time  differentiation 

X2{t,z)  =  §f  (z,i) 

partial  differentiation 

Table  6.1:  Standard  Notation  for  Lagrangian  Mechanics 


In  the  sequel,  we  refer  to  the  kinetic  energy  expression  in  terms  of  three  components; 


-frigid  body 

=  \u£lbU>b, 

(6.5) 

•^flex! 

=  5  /  l|ftw(*)  +  Vt{z)\\2pA  dz , 

Jo 

(6.6) 

^flex2 

=  k  jo[ub  +  P(t{z)]TI[u*  +  P£t{z)\p  dz. 

(6.7) 

6.1.1  FEM:  Collocation  by  Splines 


A  model  based  on  the  Finite  Element  Method  (FEM)  is  obtained  by  spatial  discretization 
via  collocation  by  splines.  The  use  of  splines  for  such  purposes  is  described  in  [BBKA88]. 
Given  the  geometric  boundary  conditions  described  above  the  discrete  approximations  of  the 
scalar  valued  functions,  Tjx(z),T]2(z),6(z),<f>(z),  are  decoupled  and  reduce  to  approximation 
of  a  single  scalar  function,  say  7 (z),  on  the  interval  z  e  [0,  l]  with  boundary  condition, 
7(0,<)  =  0.  Dividing  the  interval  [0,f]  into  N  uniform  subintervals  and  using  first-order 
B-splines  we  obtain  the  approximation: 


7(M)  «  £7i(05i-i(*)  (6-8) 


with  the  boundary  condition, 

■£Ut)BU0)  =  0. 

i= 0 

Then  following  the  reduction  procedure  described  in  [BBKA88,  §4.4.2]  we  obtain  the  FEM 
approximations  with  $(z)  as  given  in  [BBKA88,  Eqn.  (4.47)], 


ifi(M)  w  $r(z)»fr(4)  (6.9) 

rj2(z,t)  »  $T(z)fj2(t)  (6.10) 

0(z,t)  «  $r(z)0(4)  (6.11) 

<fi(z,t)  «  $T{z)4>(t).  (6-12) 
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6.1.2  Reduction  of  the  Kinetic  Energy  Function 

Using  the  FEM  approximations  of  order  N  we  obtain  the  spatial  discretization  of  the  kinetic 
energy  terms, 


where 


^flexi  ~  2  +  2wj,  fj  +  f)  IrmV} 

^flex2  ~  \{ubJ»uWb  +  2wl  -f  (  Jtil} 


l-UU)  (  )  - 


vZ  Nnr)Tj2  +  <r  -fjf  NnnT}2  ~vINn 

~ViNvvVi  ViNmV\+* 7  -V2Nn 

~Vi  Nv  ~vl Nn  Vi Nmfj2  +  fjf  N^Vi  _ 


(6.13) 

(6.14) 


(6.15) 


Nnv  =  f  pA$(z)$T(z)  dz, 
Jo 

Nn  =  f  pAz$T(z)  dz, 

[l  .  2  pAl3 
a  =  Jo  pAz  dz  =  — , 


■fuitj  - 


Nt 

1  n 


-Vi iNov  rjfN, 


*ro  J 


N, 


nv 


0 


0  N, 


Ju> u  =  f  pl{z)  dz, 
Jo 


Ju(  = 


NJ 


and 


Jii  = 


0 

0  nh 

The  following  terms  are  used  in  the  above  expressions, 


0 

Vi 

0 

Ne9 


JVf 

0 

0 


(6.16) 

(6.17) 

(6.18) 

(6.19) 

(6.20) 
(6.21) 

(6.22) 

(6.23) 


N*  =  f  pItz$T(z)  dz, 

Jo 

(6.24) 

Nj  =  f  plyyQ7^)  dZ, 

Jo 

(6.25) 

Nl  =  jT  plmx$ T{z)  dz, 

(6.26) 
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=  [' pIxx$(z)$T(z)  dz, 

Jo 

Nee  =  ftpIyy^(z)^T(z)dz. 

Jo 

Note  that  the  total  kinetic  energy  can  be  written  in  the  standard  form, 


(<jJb  ^ 


where 


mb,n,0 


lb  -f"  "h  J(t)u>  Itiiij (f?)  Ju>£ 

IIM  lor,  0 

0  J(i 


(6.27) 

(6.28) 


(6.29) 


(6.30) 


6.1.3  Reduction  of  the  Potential  Energy  and  Dissipation  Functions 

As  above,  the  potential  energy  and  dissipation  functions  are  reduced  by  substitution  of  the 
FEM  approximations.  For  the  potential  energy  we  obtain, 

1/  =  \{9TKe~0  +  <f>r K^4>  +  Vi  Kmfji  +  +  2 F Kemm  +  K^rj,},  (6.31) 


where 

Ke  =  fl{EIvyiz{z)^{z)  +  K1GA^(z)^T(z)}dz 
Jo 

K+  =  ft{EIum*z(z)**{z)  +  KtGA*(z)*T(z)}dz 
Jo 

Km  =  [\KlGA$z(z)$T,(z)}dz 
Jo 

Km  =  A«2  GA$z(z)*Tz(z)}dz 

Jo 

K«m  =  J{KiGA$(z)$*(z)}dz 

K+m  =  ~7  f  {«2 GA$(z)$*(z)}dz. 

Jo 

Similarly,  the  dissipation  function  reduces  to  the  form, 


R  =  \{fj\T Brnfji  +  fj2  Bmfj2  +  0  Be9  +  4>  B ^  +  29  B9m^ i  +  2 4>  B^rj i2}, 


_•  T 


•-.T 


where  given  the  expressions, 

Ro  :=  [  §(z)$r(z)  dz, 

Jo 

K  :=  fC*,{z)9Tiz)dz, 
Jo 

K  :=  f'i,(z)*T,(z)iz, 

Jo 


(6.32) 

(6.33) 

(6.34) 

(6.35) 

(6.36) 

(6.37) 

(6.38) 

(6.39) 

(6.40) 

(6.41) 
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we  find  it  convenient  to  express  the  individual  matrix  coefficients  in  R  (arising  from  the 
FEM  approximation)  as: 

Bm  =  Cii Ro  +  <«Ro  fori  =1,2,  (6.42) 

Be  =  (Ci  2  +  Ci3)Ro  +  Ci4-ffo>  (6.43) 

B<t>  =  (C22  +  (23)Ro  +  C24-Ro>  (6.44) 

Bme  =  -CiaR'o,  (6-45) 

B =  -C23K  (6.46) 

Finally,  we  summarize  the  model  stiffness  and  damping  using  the  expressions, 


Kr,  I  0 

0 


B  =  Bvi  0 

”  0  B„ 


.  = 


)  Bqt  — 


K»m  0 
0 


0  B i 


Kf  = 


r  B(  = 


Kg  0 
0  K+  ’ 

Bg  0 

0  5*  • 


0.1.4  Lagrange’s  Equations 

Finally,  we  compute  the  equations  of  motion  using  Lagrange’s  equations  in  the  form: 

q  =  t(q)p  (6.47) 


SV  8R  n 


(6.48) 


with  generalized  coordinates  q  =  and  velocities,  p  =  7,  £}• 


Equations  of  Motion  for  Multiaxis  Slewing:  Applying  Lagrange’s  equation  results 
(after  some  simplification)  in  the  equations  of  motion  with  kinematics  of  the  system  attitude 
expressed  in  terms  of  the  Gibbs  vector,  7; 

7  =  r(7)u%,  (6.49) 

and  kinetics, 

[ffc  "4"  Iwu  "4  Ju>u]^b  “4"  ~4  J uif. C  ”4  4“  Iu> u>  ”4  “4  “4  "4  Hlj  =  Tf,  (6.50) 

+  Irm*i  -  \HTVb  +  Rf}  +  Kr,fj  +  Kr,tl  +  Bj)  +  B^l  =  0  (6.51) 

b  +  Jtd  +  Kd  +  KvtV  +  B((  +  B*jj  =  0  (6.52) 


where 


-uinl  Nvv  -  ^Nv  2fil 

-Vinf  -Wrfjf  -  2u>3 N% 

-viNf  +  u32 \fiNm  +  fxNm  -w2Nf  -4 u,32rj\ +  fj*Nn 


R(*b) 


0  U)3N„ 


(6.53) 


(6.54) 
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6.2  Simulation  Model  for  SBL  Slewing  and  Pointing  Control  Experiments 

A  computer  simulation  model  was  developed  and  coded  based  on  the  above  FEM  model  for 
attitude  slewing  of  a  flexible  space  craft  (6.49)-(6.52).  The  model  parameters  and  geometry 
was  chosen  to  represent  the  dynamics  of  a  SBL  system  beam  expander  with  structural 
flexure  arising  from  the  flexure  response  of  the  metering  truss.  A  goal  of  the  second  year 
effort  was  to  develop  concepts  for  experimental  testing  of  nonlinear  control  concepts  for 
flexible  space  structures.  A  suitable  facility  for  testing  large  angle,  multiaxis  slewing  and 
precision  pointing  is  currently  under  construction  at  the  Astronautics  Laboratory.  The 
ASTREX  facility  with  include  3  axis  slewing  and  the  initial  test  article  is  a  scaled  model 
of  a  SBL  beam  expander.  We  have  adapted  our  nonlinear  slewing  control  simulation  model 
to  represent  available  parameters  and  geometry  of  the  ASTREX  test  article.  Details  of  the 
modeling  assumptions  and  geometry  are  given  below. 

6.2.1  Simulation  model  geometry  assumptions  and  parameters 

The  objective  of  the  simulation  model  design  was  to  generate  a  simplified  model  that  would 
duplicate,  as  closely  as  possible,  the  geometry  and  dynamics  of  the  ASTREX  test  article. 
Information  describing  the  test  article  was  sparse,  therefore,  the  following  assumptions  were 
required  to  develop  a  working  simulation  model: 

1.  The  basic  configuration  of  the  model  consisted  of  a  rigid  body,  representing  the  main 
mirror  assembly  and  associated  structure,  and  an  attached  tripod  structure,  represent¬ 
ing  the  flexible  metering  truss  with  secondary  mirror,  see  Figure  6.1. 

LOCATION  OF  ' 


Figure  6.1:  Simulation  model. 
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2.  The  metering  truss  geometry  consisted  of  a  tripod  of  tubes  located  at  the  corners  of 
an  equilateral  triangular  base.  The  generic  model  discussed  above  was  readily  adapted 
to  represent  this  geometry  by  including  spatial  variation  of  the  area  moments  Ixx(z) 
and  Iyy(z)  based  on  this  geometry.  Details  of  the  computation  of  the  area  moments 
and  the  resulting  FEM  model  approximations  are  given  in  Appendix  A. 

3.  The  dimensions  of  the  metering  truss  model  were  as  follows: 

Length,  1  -  5.0  m 

Base  width,  b  -  3.8  m 

Tube  inner  dia.  -  8.0  cm 

Tube  outer  dia.  -  12.7  cm 

Mass  of  secondary  mirror  -  190  kg  (5  %  of  total  sys.  mass) 

4.  The  overall  system  characteristics  were  as  follows: 

Mass  moments  of  inertia  Ixx  =  16640  kg  m2 

Iyy  =  16590  kg  m2 

I„  =  11660  kg  m2 

System  mass  m  =  3810  kg 

First  modal  freq.  of  metering  truss  =  12.2  Hz 

5.  The  dynamics  of  the  secondary  mirror  were  neglected  and  it  was  modeled  as  a  point 
mass  at  the  end  of  the  truss. 

6.  The  modulus  of  elasticity  of  the  truss  material  was  chosen  so  that  its  first  modal 
frequency  matched  that  of  the  ASTREX  test  article. 

0.3  Design  Slewing  and  pointing  control  with  passive  vibration  sup¬ 

pression 

The  design  methods  of  PLF  have  been  discussed  with  application  to  slewing  and  pointing  of 
a  multibody  system.  In  this  example  we  consider  the  significance  of  input-output  structure 
of  the  PLF  control  for  slewing  and  pointing  of  the  system  primary  body  axis  for  use  of  pas¬ 
sive  vibration  components.  The  essential  idea  is  that  the  use  of  passive  vibration  absorbers 
requires  tuning  for  predictable  structural  resonances.  Available  technology  for  Control  Struc¬ 
ture  Interaction  (CSI)  includes  devices  which  interact  with  the  structure  through  inertially 
generated  forces  and  torques  to  provide  vibration  absorbtion.  Such  devices  include  Proof 
Mass  Actuators  (PMA’s). 

Miller  and  Crawley  describe  and  compare  several  alternatives  for  tuning  stiffness  and 
damping  properties  of  PMA’s  for  application  as  passive  vibration  damping  of  simple  flexible 
structure  models.  Their  analysis  extends  the  classical  analysis  and  highlights  the  highly 
tuned,  narrow  band  behavior  of  such  devices  used  either  as  passive  vibration  absorbers  or 
as  active  control  actuators.  Thus  damping  effectiveness  will  depend  critically  on  assumed 
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vibration  modes.  For  nonlinear  dynamics  of  multibody  systems  such  passive  tuning  is  prob¬ 
lematical.  In  the  following  we  consider  the  importance  of  PLF  and  the  tuning  of  PMA’s  as 
passive  dampers  for  application  to  the  problem  of  rapid,  large  angle  slewing  and  precision 
pointing. 

The  slewing  and  pointing  control  design  based  on  PLF  compensation  focuses  attention  on 
stability  enhancement  of  the  decoupled  or  zero  dynamics.  Vibration  damping  should  provide 
enhanced  stabilization  in  the  sense  of  transient  response  of  the  structural  elastic  responses 
which  have  been  decoupled  from  the  primary  regulated  variables  by  the  PLF  compensation. 
Thus  PMA  tuning  should  be  optimized  for  the  system  zero  dynamics.  A  simple  example 
serves  to  illustrate  the  point.  Consider  the  case  of  slewing  and  pointing  of  the  inertial 
attitude  of  the  system  primary  body  using  external  torques  applied  to  the  primary  body. 
In  this  case  the  dynamic  response  of  the  zero  dynamics  can  be  described  by  the  (linear) 
cantilevered  vibration  modes. 

A  common  problem  in  design  of  structural  vibration  control  using  low  authority  control 
actuators  is  to  decouple  the  rigid  body  modes  from  the  vibration  control  loops  so  that  low 
authority  structural  actuators  are  not  saturated  by  interaction  with  rigid  body  dynamics. 
This  is  typically  accomplished  by  processing  paired  control  loops  from  symmetrically  placed, 
collocated  (torque)  actuators  and  sensors  [Jos89].  This  approach  may  be  sensitive  to  system 
symmetries.  An  advantage  of  the  PLF  compensation  scheme  for  primary  body  slewing  as 
described  above  is  that  the  rigid  body  motion  is  decoupled  from  the  multibody  alignment 
control. 

System  Model  for  Passive  Damping  by  PMA  The  simplest  technology  for  vibration 
absorbtion  is  based  on  vibration  damping  by  inertial  interactions.  The  Lagrangian  approach 
to  system  modeling  indicates  how  the  former  system  model  can  be  modified  to  include  the 
dynamic  interaction  between  the  PMA’s  and  the  structure. 

Consider  the  model  illustrated  in  Figure  6.2  where  two  PMA’s  with  single,  linear  degree 
of  freedom  are  located  at  the  position  z  =  t  on  the  elastic  appendage.  The  pair  of  PMA’s 
are  assumed  to  have  axes  of  motion  aligned  with  the  local  coordinate  frame.  The  motion  of 
the  moving  mass  for  PMA  aligned  with  the  ®-axis  (resp.  y-axis)  is  given  by  the  coordinate 
x„  (resp.  y„)  defined  with  respect  to  a  zero  energy  at  rest  configuration. 

The  modified  system  Kinetic  Energy  has  the  form, 

^  =  frigid  body  -^flex  +  Tpma- 
The  inertial  velocity  of  the  appendage  at  a  point  z  is: 

l  \ 

v(z)  :=  t?(z)  +  tj,(z) 

W(z)  +  z) 

Neglecting  axial  deformation,  773  «  0,  we  obtain  the  additional  kinetic  energy  in  the  form, 

Tpma  =  \rna{xa  -  tq(z))2  4-  Jma(y0  -  v2{z))2  (6.55) 
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Figure  6.2:  SBL  model  with  passive  damping  by  PMA’s. 

=  \rna[xa  -  tji(z)  -  u2z  +  u9rj2(z)]2 

+kma[ya  ~  rii(z)  +  u/i z  -  u9t)i(z))2. 

Assuming  the  devices  are  identical  and  decoupled  we  take  the  system  potential  energy  in  the 
form, 

V  =  VfLex  +  vpma, 

where, 

VpMA  =  \KXa  -I'  \KVa2 ■ 

Applying  Lagrange’s  equations  to  the  PMA  terms, 

It  (^r)  ~  ^ Tpua  ~  VpMi)  =  Q " 

with  respect  to  the  expanded  set  of  reduced  generalized  coordinates  obtained  by  retaining 
N  finite  elements, 

€  {lb,Vi,rj2,0,<t>,xa,ya}, 

one  obtains  the  modified  equations  of  motion  in  the  standard  form: 

q  =  f  {q)p  (6.56) 

M,y,p  +  B.y,p  +  Ktylq  =  G,y,n  (6.57) 

where 

_  _  IT  -T]T 

9  =  I7n  *?.?»*  J  » 

r  T  -T  jp  ,  ' 

v  =  K » v  ,  £  .  *  ]  i 

and  *  =  [®a,  ya]T.  Denoting  terms  arising  from  the  interaction  with  the  PMA’s  by  superscript 
‘a’  we  can  write  the  modified  N ,  x  N,  system  matrices  (with  Nt  =  2 N  +  5)  in  the  form, 


M,y.  — 


h  +  IUU  +  Jww  +  Af“w 

lu >TJ  4" 

Jut 

AT® 

WX 

IL,  +  [Mi/ 

■Ljtj  + 

0 

M“ 

lit 

0 

Jtf 

0 

(ACf 

0 

Ma 
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B.y s  = 


£lb{h  +  Iww  +  Jwu>  + 

fife  l/ 

b, ' 

IT  +  Ba 

-  urn  '  wn 

i2  +  5,,  -f- 

0 

■s. 

0 

0 

TDa 

Dpma  . 

K,y,= 


'  0 

' 

Kr, 

K,( 

0 

*1 

Kt 

[  0 

Kpma 

,  and  Gaya  — 


h 

0 

0 

0 


Letting  both  PMA's  be  located  at  the  end  of  the  appendage;  z  =  l,  then  the  new 
take  the  form, 


m„£2 


0 


Ma 


M“  = 

U?T) 


m:x  =  - 


—matrf\{t) 

0  ma£2  —ma£7]2(l) 

[  -ma£t]x(£)  -mji]2(e)  ma(i}\(£)  +  tjI(£))  J 


(3  x  3) 


0---0  0 
0  •  •  •  0  ma£ 

L  0  •  -  -  0  -maT]2(£) 

0 

ma£ 

-mati2{£) 


0-**0  —ma£ 

0--.0  0 
0*-'0  -mar]i(£)  J 


(3  x  2 N) 


-ma£ 

0 

-rnarji(£) 


(3x2) 


M“  = 

TJX 


'  0 

* 

• 

0 

0 

— ma 

0 

0 

0 

-ma 

(2 N  x  2) 


=  diag{0,...,0,ma,0,...,0,mo}  (2N  x  2N) 


B* 

IJV 


Ba 


12  — 


0 

Bh] 

—  Ra 

aX2 

0 

(2 N  x  2 N) 


0 

0 

• 

0 

0  •••  0 

-2mau>3 
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PMA  Tuning  for  Slewing  and  Pointing  Control  Such  devices  provide  dissipation  by 
coupling  kinetic  vibration  energy  into  the  moving  mass  of  the  PMA  where  it  can  be  dissipated 
by  localized  damping.  Since  the  effectiveness  of  the  energy  coupling  depends  on  inertial 
interactions  the  mass  of  the  device  is  a  limiting  factor  in  its  effectiveness.  Nevertheless,  the 
effectiveness  of  the  passive  damping  depends  on  tuning  the  PMA  stiffness  and  damping  to 
the  critical  vibration  mode. 

Miller  and  Crawley  review  classical  and  modern  approaches  to  tuning  PMA’s  as  vibration 
absorbers.  They  discuss  tuning  with  respect  to  the  simple  linear  model  of  a  PMA  interacting 
with  a  single  degree-of-freedom  system  as  shown  in  Figure  6.3  where  mx  and  kx  are  the 
structure  mass  and  stiffness  respectively,  and  where  m2,  fc2,  and  c  are  the  absorber  mass, 
spring  stiffness,  and  damper  constant  respectively.  We  focus  attention  on  tuning  for  optimal 
transient  response  characteristics  as  follows.  Given  the  non-dimensional  parameters: 


6 

0 


\ 


{k2 /m2) 

(ki/mi) 


absorber  -  structure  frequency  ratio 


m2  ,  , 

—  absorber  -  structure  mass  ratio 
mi 

non-dimensional  damping 


2m2\JkiJrn\ 


choose  the  PMA  stiffness  ka  and  damping  c  for  a  fixed  mass  ratio  such  that  the  maximum 
modal  time  constant  is  minimized.  The  tradeoff  is  illustrated  in  Figure  6.4  which  displays  a 
root  locus  as  fi  is  varied  (for  a  mass  ration  /3  =  0.02).  The  frequency  ratio,  6 ,  for  optimal 
transient  response  is  the  value  for  which  the  root  loci  intersect  at  point  A,  and  is  given  by, 

Str  =  Y+jjj'  (6-67) 

Application  of  this  approach  to  flexible  structures  is  of  course  limited  by  the  interaction 
between  multiple  vibration  modes. 

For  the  problem  of  slewing  and  pointing  as  discussed  in  the  previous  reports  we  let  the 
primary  outputs  be  the  inertial  attitude  parameters,  -y*,,  of  the  primary  body  frame  and 
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c 


Figure  6.3:  Model  used  by  Crawley  and  Miller  for  tuning  passive  absorber  to  a  one  DOF 
system 


-O.iO  -0.05  0.00 

REAL  PART 


Figure  6.4:  Root  locus  of  poles  for  Miller  and  Crawley’s  one  DOF  system  as  a  function  of 
increasing  damper  value 
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Mode 

Frequency 

1 

12.07  Hz 

2 

12.07 

3 

70.73 

4 

70.73 

5 

103.19 

6 

103.19 

7 

271.33 

8 

271.33 

Table  6.2:  Truss  modal  frequencies  for  cantilevered  boundary  conditions 

apply  external  torques,  77,,  to  the  primary  body.  In  this  case  the  decoupled  dynamics  are 
given  by  the  linear  dynamics  of  the  metering  truss  with  cantilevered  boundary  conditions 
replacing  the  interaction  with  the  rigid  body.  (To  see  this  simply  constrain  the  primary 
regulated  outputs,  y  =  7^  and  their  rates  y  to  the  origin  and  substitute  7^  =  0,  u>(,  =  0  in  the 
equations  of  motion  to  obtain  the  zero  dynamics.)  To  tune  the  PMA’s  for  optimal  transient 
response  of  the  first  pair  flexible  modes  of  the  zero  dynamics  we  apply  the  strategy  of  Miller 
and  Crawley  to  the  cantilevered  appendage  model. 

Utilizing  the  results  of  Miller  and  Crawley,  we  chose  an  absorber  mass  to  truss  mass 
ratio  of  .05  for  each  device  and  performed  a  root  locus  analysis  with  Str  given  in  (6.67). 
The  root  locus  appears  in  Figure  6.5  showing  the  poles  of  the  passive  absorber  and  lowest 
modal  frequency  of  the  truss  for  different  values  of  absorber  damping.  Maximum  structural 
damping  occurred  at  the  intersection  of  the  root  trajectories,  with  optimal  values  of  absorber 
damping  and  spring  stiffness  of  370  kg/s  and  3.9  xlO *N/m.  As  expected,  the  optimal 
frequency  ratio  suggested  in  (6.67)  offered  a  crude  approximation  and  further  refinement 
was  necessary  to  obtain  transient  response  tuning  even  for  the  simple  finite  element  model 
used  in  our  simulations.  The  resulting  tuned  values  of  PMA  parameters  were  then  used  in 
simulations  of  PLF  slewing  and  pointing  (refer  to  model  MDDA). 


Summary  of  Simulation  Results  Computer  simulations  of  the  finite  element  model 
HODA  were  performed  comparing  the  effect  the  passive  PMA  absorber  on  slew  times- and 
torque  requirements.  A  standard  three  axis  slew  maneuver  was  defined  by  initial  offset  from 
the  desired  inertial  attitude  given  by  the  Gibbs  vector, 

76  =  [0.08,0.09,0.036] 


or  in  Euler  angles, 


=  [5.°,10.o,10°] 


and  results  are  summarized  in  tables  6.3  and  6.4.  Slew  time  is  computed  as  the  elapsed  time 
for  the  absolute  value  of  each  Gibbs  component  of  the  rigid  body  attitude,  7s,  to  become 
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Figure  6.5:  Root  locus  showing  interaction  between  absorber  poles  and  lowest  frequency 
truss  poles  as  a  function  of  absorber  damping 

less  than  0.0002.  Peak  torque  requirements  are  computed  relative  to  a  fixed  torque  level, 
r0,  equivalent  to  angular  acceleration  of  10  deg/sec/sec  for  an  equivalent  rigid  body  with 
inertia  moment  equal  to  the  total  undeformed  ASTREX  system  inertia.  The  peak  torque  is 
computed  as 

T-peafc  =  max{f  >  0  :  - — ■  }■ 

lro| 

The  tables  also  include  performance  tradeoffs  for  implementation  of  PLF  for  various 
reduced  order  models  with  r  retained  degrees  of  freedom  in  the  SBL  system  model.  Results 
indicate  that  slew  times  were  improved  with  PM  A  tuned  for  the  dominant  elastic  mode 
of  the  cantilevered  truss  for  various  reduced  order  PLF  implementations.  This  confirms 
the  theoretical  predictions  of  the  PLF  compensation  and  stability  of  the  zero  dynamics. 
However,  peak  torque  requirements  remained  relatively  unchanged.  (Note  that  in  Tables  6.3 
and  6.4,  MODA  r=6  and  r=7,  correspond  to  M0DB,  r=4  and  r=5,  since  M0DA  r=4  and  r=5  are 
the  modes  associated  with  the  absorber  dynamics  which  are  not  included  in  MODB.).  This 
follows  by  observing  the  additional  mass  of  the  PMA  offset  their  benefit  in  stabilization 
with  reference  to  slewing  torque  requirements.  Figures  6.6  and  6.7  show  time  responses  for 
slewing  simulations  comparing  SIHA5,  a  reduced  order  model  with  absorber  which  retains 
the  absorber  modes  and  first  two  truss  modes,  and  SIMC2,  a  reduced  order  model  without 
absorber  which  retains  first  two  truss  modes.  The  dotted  traces  correspond  to  simulation 
SIMAS  and  illustrate  the  dramatic  increase  in  truss  damping  generated  by  the  tuned  absorber 
as  compared  with  the  solid  traces  of  simulation  SIMC2.  (See  also  Section  6.3  for  summary  of 
simulation  results.) 
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Peak  Torque  ATpk(r)/Tpk(Z ) 

Relative  to  rigid  body  case 

simulation 

PLF  Model  DOF 

model  r  = 

H 

4 

5 

6 

7 

13 

0 

0.01 

0.48 

0.89 

111 

KSlfli 

0 

0.86 

EH 

Table  6.3:  Relative  Peak  Torque  Increase  with  PLF  Compensation  DOF 


Slew  Times  (sec)  to  |-y»|  <  .0002 


PLF  Model  DOF 

r  = 

3 

4 

5 

6 

7 

13 

Model  MODA 

Axis  1 

Wjpfpl 

1.472 

1.471 

1.419 

2 

Baa 

nil 

1.627 

1.483 

1.476 

3 

0.940 

BH 

m 

0.940 

Mode 

K0DB 

Axis  1 

1.412 

2 

1.638 

1.638 

1.482 

1.469 

3 

0.941 

Table  6.4:  Slewing  Times  for  Standard  3-axis  Maneuver 
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Figure  6.6:  Time  responses  comparing  and  rj2(t)  for  the  models  with  and  without  tuned 
absorber 


Figure  6.7:  Time  responses  comparing  0(t)  and  (f>(( )  for  the  models  with  and  without  tuned 
absorber 
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0.4  Design  Example:  LOS  Slewing  and  Multibody  Alignment 

The  application  of  slewing  and  pointing  control  to  SBL  as  well  as  other  space-based  optical 
systems  will  require  the  coordination  of  multiple  actuators  for  system  slewing  and  multi¬ 
body  alignment.  The  effective  optical  LOS  for  a  multiple  mirror  system  will  involve  relative 
motion  of  the  mirrors  as  well  as  the  primary  body  inertial  attitude.  Detailed  analysis  of  the 
optical  train  of  a  candidate  SBL  system  is  beyond  the  scope  of  this  study.  We  will  consider 
a  simplified  problem  which  highlights  a  significant  aspect  of  the  application  of  PLF  compen 
sation  for  multibody  alignment  when  the  available  controls  and/or  the  choice  of  regulated 
outputs  leads  to  system  zero  dynamics  which  are  unstable  in  some  sense.  Of  course,  the 
sense  in  which  instability  of  the  zero  dynamics  may  limit  application  of  PLF  compensation 
depends  on  the  application.  The  available  theory  suggests  that  PLF  compensation  is  not 
feasible  for  nonminimum  phase  systems. 

Clearly,  the  limitation  of  applying  PLF  compensation  to  a  nonminimum  phase  system 
arises  from  the  fact  that  PLF  compensation  seeks  to  implement  an  online  system  inverse 
which  is  by  definition  unstable  if  the  system  model  is  nonminimum  phase.  For  the  class  of 
Lagrangian  systems  arising  from  control  of  flexible  structures  the  PLF  compensation  includes 
an  Inverse  Force  Model.  In  this  section  we  illustrate  some  tradeoffs  in  the  implementation 
of  an  online  IFM  for  a  nonminimum  phase  system. 

System  Model  for  Multibody  Alignment  Control:  The  method  of  Lagrange’s  equa¬ 
tions  offers  a  systematic  approach  to  develop  equations  of  motion  for  systems  with  external 
forces  described  by  a  virtual  work  function.  We  consider  the  two  body  system  described  pre¬ 
viously  where  we  have  assumed  the  multibody  configuration  rotates  about  the  point  z  =  0 
at  which  the  appendage  is  attached  to  the  primary  body.  We  now  assume  that  the  following 
external  forces  are  applied  to  the  two-body  system: 

1.  torques,  Tfc,  applied  to  the  primary  body, 

2.  distributed  torque  function  applied  to  the  flexible  appendage,  r(z)  =  Y,kTkb{z  —  Ck), 

3.  distributed  force  function  applied  to  the  flexible  appendage,  f(z)  =  J2j  fj$(z  ~  C i )• 

The  virtual  work  function  can  then  be  derived  as  follows.  The  power  flow  into  the  system 
is  given  as, 

P  =  rigid  body  torque 

+  /orT(z)[o;j,  +  |(z)J  dz  distributed  torque 
+  Jo  fT(z)lub  xz  +  r}(z)]  dz  distributed  force 

where  z  =  [0, 0,  z)T.  Note  that 

’  0  —  z  0 

u x  2  =  -2  x  wj  =  -  z  0  0 

0  0  0 
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and  since  u>b  =  r(£j,)£i  we  can  write 


P  =  T?nM(>  +  j‘{TT(*mt)(l,  +  TT(z)i(z)}dz  (6.68) 

+ J‘{-fT(z)zr((t)(t  +  fT(z)i(z)}  iz. 

The  virtual  work  in  (6.4)  is  replaced  with, 

8W  =  |r,rr(6)  +  foTT{z)T(ib)dz-  Jl  fT(z)ZT(tb)dz^  6(b  (6.69) 

+  {jfE^(*-Ck)  H  dz J 

+  |^  X]  fj S(z  -  GO  Srl  dz\  . 

Defining  the  vector  of  external  forces, 

/.  =  [fu-  (6.70) 

and  torques, 

f,  =  [r1,...,rfe,...]r,  (6.71) 

and  applying  Lagrange’s  equations  gives  the  system  equations  of  motion  in  the  form, 

[lb  +  hu  +  Jww]d>b  +  Iutf  +  +  /w  +  Juu]vb  +  Iwvfj  +  -f  Hfj 

=  Tb  +  'Erk  +  '£Zjfj  (6.72) 

fc  i 

+  Innf}  -  -I-  Rij  +  Kvfj  +  +  Brjf}  +  Br,^,  =  Gj,  (6.73) 

J^b  +  +  K*(fj  +  Bt(  +  B^fj  =  G(T,  (6.74) 

where  H,  and  Jl  are  given  by  (6.53)-(6.54)  and 


Then  letting 


0 

-01 

II 

NT 

0 

0 

. 

0 

0 

■  1  0 

1  0 

- 

- 

0  1 

. . . 

0  1 

. . . 

.  0  0 

0  0 

z 

=  [Zu 

Zi\  ■ 

•  >  Zjy 

...] 

we  can  express  the  equations  of  motion  in  the  standard  form  (6.56)-(6.57)  with 


G  = 


‘hi  Z 
0  0  Gr, 

.  0  C?<  0 


(6.75) 
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A  Control  Configuration  for  System  Slewing  and  LOS  Pointing  Slewing  and  point¬ 
ing  control  of  the  optical  LOS  of  a  typical  beam  expander  for  a  SBL  system  will  require 
system  slewing  and  alignment  of  the  primary  and  secondary  mirrors  of  the  optical  system. 
Such  requirements  are  typical  of  various  large  aperture  optical  systems.  With  reference  to 
our  two-body  SBL  model  the  primary  body  inertial  attitude — given  by,  in  Euler  angles, 
(or  7 in  Gibbs  parameters) — represents  the  system  pointing  of  the  primary  mirror.  Op¬ 
tical  misalignments  contributing  to  LOS  pointing  errors  can  be  described  by  misalignment 
variables ; 


2/4  =  9(t)  +  slVl{£)  (6.76) 

t/5  =  <t>{i)  -  slV2{t)  (6.77) 

where  $i  is  a  sensitivity  parameter  which  depends  on  the  optical  components.  Then  slewing 
of  the  system  LOS  can  be  viewed  as  a  problem  of  slewing  the  primary  body  attitude  while 
regulating  the  misalignment  variables  to  zero. 

Nonminimum  Phase  Characteristics  and  Actuation  We  assume  that  control  is  im¬ 
plemented  by  three  independent  external  torques  applied  to  the  primary  system  body  to¬ 
gether  with  two  independent  control  forces  applied  at  the  secondary  body.  To  implement 
PLF  control  for  the  resulting  5-input /5-output  problem  we  first  identify  stability  characteris¬ 
tics  of  the  zero  dynamics.  Since  our  concern  is  with  regulation  of  the  primary  system  outputs 
consisting  of  the  system  outputs  &  and  the  misalignment  variables  and  stability  of  the  elastic 
structure  with  respect  to  its  undeformed  state  we  focus  attention  on  the  stability  properties 
of  the  linear  perturbation  model  for  the  flexible  structure  as  &  =  0  and  y4  =  yB  =  0.  In 
this  case,  the  local  stability  of  the  zero  dynamics  in  a  neighborhood  of  the  undeformed  state 
can  be  determined  from  the  zeros  of  the  transfer  function  describing  dynamics  of  the  linear 
perturbation  model. 

The  sensitivity  parameter  which  defines  the  misalignments  contributes  to  the  stability 
properties  as  shown  Figure  6.8.  The  zero  locus  shows  a  nonminimum  phase  condition  occurs 
for  si  small  enough.  The  parameter  si  can  be  seen  as  a  sensitivity  factor  that  proportions 
the  importance  of  the  linear  deformations  relative  to  the  angular  deformations  in  the  outputs 
and  indicates  that  it  is  difficult  to  control  angles  using  forces  instead  of  torques. 

Practical  Aspects  of  Implementing  PLF  Compensation  for  Nonminimum  Phase 
Systems  The  above  simple  model  indicates  how  parametric  sensitivities  can  change  the 
nature  of  the  nonlinear  control  problem  considered.  Available  theory  for  PLF  compensation 
is  predicated  on  the  assumption  that  the  system  is  nonminimum  phase.  At  this  point  in  our 
study  we  were  interested  in  considering  practical  situations  which  might  be  encountered  in 
experiments  which  violated  the  assumptions  which  might  limit  application  of  PLF  control. 
We  wanted  to  investigate  several  questions:  To  what  extent  model  order  reduction  in  the 
design  of  PLF  compensation  is  effected  by  nonminimum  phase  characteristics?  Systems  with 
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Poles  zeros  -  linear  coab  of  theta  etal,  phi  eta: 


Figure  6.8:  Locus  of  zeros  for  model  MODE  as  Si  is  increased 

very  fast  unstable  zero  dynamics  modes  will  be  characterized  by  linear  models  with  right  half 
plane  zeros  which  are  far  from  the  imaginary  axis.  Can  model  reduction  schemes  be  adapted 
to  remove  such  characteristics  from  the  nominal  design  model?  If  PLF  compensation  is  then 
designed  for  the  reduced  and  minimum  phase  model,  can  it  be  applied  reliably  to  the  full 
order  system? 

For  implementation  and  robustness  evaluation  of  PLF  control  we  chose  a  value  of 
which  gives  nonminimum  phase  behavior.  We  first  applied  modal  truncation  to  the  open 
loop  system  model  retaining  various  order  models  to  determine  if  a  nonminimum  phase  model 
could  be  found.  The  poles  and  zeros  of  various  reduced  order  models  using  modal  truncation 
appear  in  Tables  6.5  and  6.6.  The  number  of  retained  dynamic  degrees  of  freedom  is  indicated 
as  r.  For  r  =  6  the  linear  perturbation  model  has  no  right  half  plane  zeros.  However,  for 
r  =  7  a  right  half  plane  zero  appears  indicating  a  nonminimum  phase  condition.  Simulation 
studies  were  performed  to  test  the  behavior  of  the  nonlinear  system  response  with  PLF 
compensation  implemented  based  on  the  reduced  order  system  model.  Cases  with  5  <  r  <  7 
revealed  stable  responses. 

0.5  Design  Example:  Slewing  Control  with  Noncollocation  of  Control  and  Pri¬ 
mary  outputs 

In  the  ASTREX  test  article  configuration,  the  rigid  body  torques  are  to  be  generated  by 
cold  gas  thrusters  located  a  distance  from  the  center  of  mass  of  the  system,  which  is  the 
location  of  the  attitude  sensors.  To  simulate  the  effects  of  this  non-colocation  configuration, 
simulations  model  MODD  and  MODF  were  designed  so  that  forces  Fx  and  Fy  at  the  end  of  the 
truss  are  to  be  used  to  control  74.  The  torque  produced  by  these  forces  was  limited  to  the 


69 


Techno-Sciences,  Inc. 


TSI-90-10-10-WB 


70 


Techno-Sciences,  Inc. 


TSI-90-10-10-WB 


maximum  torque  available  on  ASTREX  for  attitude  control.  These  maximum  torque  values 
were  determined  from  the  thruster  specifications  and  the  moment  arms  measured  relative  to 
the  x,  y,  and  z  axes.  The  thruster  locations  are  shown  in  Figure  6.9,  and  each  produces  a 
maximum  thrust  of  889.6  N.  This  configuration  yields  the  following  maximum  body  torques 
which  were  used  for  all  simulations  except  those  of  model  MODD: 

Tmaxz  -  8600  Nm 
Tmax ,  =  4964  Nm 
Tmaxr  =  4964  Nm 

The  equivalent  maximum  forces  used  for  Fx  and  Fy  in  simulations  SIMD1-SIMD5  appear 
below. 


Vmoi 


992.8  N 
1720.0  N 


A  pole-zero  plot  of  the  linear  perturbation  model  of  this  system  appears  in  Figure  6.10 
and  shows  a  stable  and  nonminimum  phase  response.  Simulations  for  various  reduced  order 
models  for  a  maneuver  of 


76  =  [0.05,0.05,0] 


(or  in  Euler  angles, 


[^>0,0]  =  [0.3°,  5.7°,  5.7°]) 


appear  in  Appendix  A.  Various  model  order  reductions  were  tested  for  implementation  of 
PLF  compensation  for  slewing  of  the  primary  body  attitude  and  regulation  of  the  misalign¬ 
ments. 


Rapid  Slewing  with  multiple  modes  of  actuation.  As  discussed  in  the  report  [Benn89] 
the  use  of  combined  continuous  and  discontinuous  modes  of  actuation  can  be  incorporated 
in  the  rapid  slewing  and  precision  pointing  problem.  Focusing  attention  on  the  physical 
limitations  of  the  ASTREX  facility  we  performed  simulation  studies  of  primary  body  slewing 
as  above  but  operating  the  applied  forces  at  the  end  of  the  appendage  in  discontinuous  mode. 
In  a  physical  system  the  applied  forces  might  be  obtained  using  reaction  jets..  A  critical 
technology  being  developed  for  the  ASTREX  facility  involves  throttleable  cold  gas  jets.  In 
our  model  we  have  assumed  the  jets  can  be  modeled  as  providing  pure  external  forces  which 
can  change  subject  to  slew  rate  limiting. 

6.0  Summary  of  Simulation  Results  for  Rapid  Slewing  and  Precision  Pointing 

A  computer  simulation  of  the  SBL  two-body  model  was  developed  in  PC-MATLAB  to  eval¬ 
uate  slewing  and  LOS  pointing  performance.  The  simulation  model  parameters  shown  in 
Table  6.7  are  taken  roughly  from  available  data  on  the  ASTREX  test  article. 
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As  can  be  seen  from  the  table,  five  simulation  models  with  three  input/output  configu¬ 
rations  were  developed.  The  purpose  of  models  MODA,  MODB,  and  MODC  was  to  examine  how 
slewing  performance  was  affected  when  passive  absorbers  were  added  to  the  end  of  the  truss. 
These  models  were  configured  with  the  control  inputs  as  rigid  body  torques,  r&,  at  each  axis 
and  measured  outputs  as  the  rigid  body  attitude,  7^,  in  Gibbs  vector  form.  Figure  6.11 
illustrates  the  input  output  configuration  for  these  models.  Also,  two  passive  absorber  de¬ 
vices,  described  in  detail  above,  are  located  at  the  end  of  the  truss  for  vibration  damping. 
The  three  models  differed  only  in  the  values  of  the  parameters  chosen  for  the  spring  and 
damper  constants.  In  Model  MODA,  the  parameters  were  tuned  for  optimal  transient  response 
as  described  above.  In  model  MODB,  the  passive  damper  was  tuned  to  twice  the  highest  MODAl 
frequency  to  put  it  in  a  de-tuned  configuration,  and  model  MODC  the  passive  damper  was 
decoupled  from  the  truss  by  making  both  the  spring  and  damper  constant  equal  to  zero. 
Tradeoff  studies  comparing  simulation  results  of  these  models  were  performed. 

The  purpose  of  model  MODD  was  to  examine  how  non-colocated  sensing  and  actuation 
affects  slewing  performance.  This  is  important  because  the  ASTREX  test  article  is  configured 
such  that  body  torques  are  generated  by  thrusters  that  are  not  located  at  the  center  of 
mass  of  the  body  which  is  the  location  of  the  attitude  sensors.  In  order  to  simulate  this 
non-colocated  condition,  model  MODD,  shown  in  Figure  6.12  was  set  up  with  control  inputs 
as  two  force  actuators,  Fx  and  Fy,  at  the  end  of  the  truss  and  a  torque  t3,  about  axis 
three  and  measured  outputs  as  the  rigid  body  attitude  ,7 &,  in  Gibbs  vector  form.  Although 
this  configuration  is  not  an  exact  representation  of  the  thruster  placement  on  ASTREX,  it 
provided  us  with  an  approximate  model  to  study  this  non-colocated  condition. 

Finally,  the  purpose  of  model  MODE  was  to  examine  slewing  performance  and  stability  of 
the  SBL  with  additional  force  actuators  added  at  the  end  of  the  truss  to  control  its  vibration. 
Model  MODE,  shown  in  Figure  6.13  was  set  up  with  the  control  inputs  as  two  force  actuators, 
Fx  and  Fy,  at  the  end  of  the  truss  and  rigid  body  torques,  -q,,  at  each  axis  and  measured 
outputs  as  the  rigid  body  attitude,  7&,  in  Gibbs  vector  form  and  two  additional  outputs 
relating  to  the  line  of  sight  error. 

Table  6.8  identifies  the  simulations  performed  for  each  model  along  with  the  correspond¬ 
ing  control  law.  Each  simulation  is  summarized  in  Appendix  A  by  plots  illustrating  responses 
of  position,  7 &,  attitude  rate,  7*,,  primary  body  torques,  q,,  any  external  forces,  Fx,Fy  ,  etc., 
any  other  outputs,  and  deflections  at  the  end  of  the  truss,  7(f),  £(f). 

0.7  Experimental  Protocol  for  Validation  of  PLF:  Mode  Locking  Experiments 

A  principal  objective  of  the  proposed  experiments  is  to  verify  the  effectiveness  of  PLF  com¬ 
pensation  by  validating  the  extent  to  which  distinctly  nonlinear  behavior,  observable  in  the 
system  principal  input-output  responses  for  attitude  slewing  and  pointing,  can  be  compen¬ 
sated.  Our  view  is  that  the  occurrence  of  mode  locking  (sometimes  called  phase  locking  or 
entrainment  [GH83])  offers  a  framework  for  observing  nonlinear  response  of  pointing  and 
tracking  control  loops  in  experiments.  Such  nonlinear  behavior  can  always  occur  when 
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Models 

Control  Laws 

|  Continuous  PLF/Slewing 

j  Discontinuous  PLF/Slewing 

r  =  3 

HI 

r  =  5 

HI 

HI 

HI 

Slew  rate 
limited  torque 

Discontinuous 

torque 

Multiaxis  Slew 

mim 

m 

SIMA3 

SIMA4 

SIMA5 

SIMA6 

SIMA7 

mm 

SIMB3 

SIMB4 

■ 

Baa 

SIMC2 

SIMC3 

MODD 

SIMD1 

1 

SIMD2 

SIMD3 

SIMD4 

SIMD5 

MODE 

SIME1 

SIME2 

Table  6.8:  Matrix  of  PLF/Slewing  Control  Implementations 


conventional  linear  attitude  control  design  is  applied  for  muitiaxis  attitude  control  imple¬ 
mentation,  but  can  be  removed  by  the  use  of  PLF  compensation.  Mode  locking  is  only  one  of 
several  types  of  nonlinear  phenomenon  which  might  be  considered  for  experimental  demon¬ 
stration.  We  propose  to  focus  experiments  on  demonstration  of  mode  locking  for  several 
reasons:  mode  locking  is  easy  to  observe  in  experiments,  is  repeatable,  and  a  distinctly  non¬ 
linear  phenomenon.  Mode  locking  is  also  of  fundamental  significance  to  achievable  accuracy 
of  precision  poking  control  laws  for  space  craft  attitude  both  with  and  without  flexible, 
multibody  interactions. 

Mode  locking  can  be  readily  demonstrated  in  attitude  slewing  and  pointing  experiments 
with  the  ASTREX  facility  and  initial  test  article  since  it  requires  relatively  small  angle  mo¬ 
tions  and  can  (in  principal)  be  demonstrated  in  various  frequency  bands.  The  experiments 
require  periodic  excitation  of  the  system  attitude  which  permits  focusing  attention  on  the 
frequency  response  and  its  dependence  on  elastic  interactions  at  higher  frequencies.  This  is 
significant  for  the  experimental  protocol  we  have  planned  in  which  we  emphasize  the  relation¬ 
ship  between  control  law  precision  and  closed  loop  bandwidth.  Multibody  interactions  and 
structural  flexure  will  complicate  the  types  of  mode  locking  possible  but  its  demonstration 
in  experiments  should  still  be  easy  to  observe. 

Moreover,  mode  locking  is  a  relevant  factor  in  control  system  performance.  It  is  well 
known  that  certain  nonlinear  affects  can  significantly  degrade  feedback  system  performance 
such  as  accuracy  of  tracking.  Mode  locking  will  contribute  to  a  degradation  of  closed  loop 
tracking.  From  the  perspective  of  nonlinear  dynamics  it  is  well  known  that  mode  locking 
is  often  a  precursor  to  chaotic  behavior  [BBJ84]  which,  in  the  case  of  closed  loop  tracking, 
would  imply  a  complete  loss  of  tracking  effectiveness. 

In  principle,  PLF  can  be  used  to  compensate  for  the  nonlinear  interactions  which  cause 
mode  locking,  as  well  as  other,  more  complicated,  dynamical  behavior.  Thus,  feedback 
linearization  should  be  viewed  as  a  means  for  introducing  nonlinear  compensation  in  the 
design  of  feedback  loops  for  achieving  robust  precision  tracking  as  well  as  executing  large 
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angle,  multi- axis  slewing. 

Mode  Locking  Experimental  Protocol.  We  next  describe  in  general  terms  the  proto¬ 
col  for  experimental  demonstration  of  mode  locking  in  the  context  of  attitude  control  and 
compensation  by  PLF. 

Mode  locking  is  a  phenomenon  associated  with  nonlinear  dynamics  involving  interactions 
of  periodic  motions  with  at  least  two  distinct  frequencies.  The  underlying  periodic  behavior 
may  be  self  induced  or  associated  with  external  excitation  as  in  active  control  loops.  The 
simplest  and  most  common  example  of  a  system  exhibiting  mode  locking  is  the  periodically 
forced  Van  der  Pol  oscillator  [Sto50,  GH83,  TS86].  In  this  case,  the  system  exhibits  a 
self  excited  oscillation  (limit  cycle).  Many  other  periodically  forced  and/or  interconnected 
oscillators  have  been  studied  using  both  analytical  and  experimental  demonstrations  of  mode 
locking  [Sto50,  GH83,  TS86,  BBJ84], 

One  experiment  would  involve  the  design  and  testing  of  multiaxis  attitude  tracking  con¬ 
trol.  The  experiment  would  involve  driving  the  attitude  control  loops  with  two  independent 
periodic  command  signals.  To  demonstrate  small  signal  tracking  response  of  the  attitude 
control  loops  one  could  use  command  signals  as, 

4>c(t)  =  A\  sina>i<,  0c(t)  =  A2  sin  w2t, 

whereas,  to  demonstrate  slewing  response  the  commands  would  be  replaced  with  their  square 
wave  equivalent.  In  either  case,  we  anticipate  mode  locking  to  be  observed  in  the  following 
way.  First,  assume  that  A2,w2  are  fixed  and  Ai,w\  are  adjustable  parameters.  As  >1!  — ►  0 
we  expect  to  observe  a  periodic  oscillation  in  the  system  attitude  response  of  frequency  u/2. 
Clearly,  if  A2  is  small  the  system  response  is  near  that  of  a  linear  system  which  can  be 
confirmed  in  the  experiment  as  follows.  Define  the  ratio,  0  :=  ^  and  choose  integers,  p,  q 
with  q  >  p,  then  the  frequencies  are  said  to  be  commensurable  if  they  satisfy  the  relation, 

pu>\  —  qu>  2  =  0 

and  fi  =  p/q  is  a  rational  number.  For  sufficiently  small  driving  amplitude  A\,  the  system 
attitude  response  will  be  periodic  with  frequency  u/2  and  q  is  a  subharmonic  number  asso¬ 
ciated  with  this  periodic  motion.  For  linear  systems,  if  f l  is  nor  rational  then  the  attitude 
motions  are  not  periodic. 

For  the  nonlinear  system,  if  Q  is  not  rational  and  Ai  is  small,  the  motion  will  not  be 
periodic  (as  in  the  linear  case).  However,  if  the  driving  amplitude  A\  is  increased  then — after 
a  transient — the  motion  will  again  become  periodic.  This,  of  course,  is  distinctly  nonlinear 
behavior.  The  periodic  motions  are  called  mode  locked  (or  simply  locked).  The  aperiodic 
motions  are  referred  to  as  drifting  behavior.  Such  behavior  is  tj'pical  of  nonlinear  systems 
and  can  be  readily  demonstrated  in  experiments.  It  is  common  to  display  mode  locking 
responses  as  regions  of  locked  and  drift  behavior  in  a  (j4,fi)  parameter  space  diagram  as 
shown  in  Figure  6.14.  Periodic  (locked)  and  aperiodic  (drift)  can  be  readily  identified  in 
measured  attitude  responses  observed  in  experiments. 


77 


Techno-Sciences,  Inc. 


TSI-90-10-10-WB 


Figure  6.14:  Typical  parameter  space  diagram  showing  regions  of  locked  (shaded)  and  drift 
behavior. 

We  expect  that  mode  locking  will  be  observable  when  a  conventional  attitude  controller 
is  used  because  the  nonlinear  kinematic  and  dynamic  couplings  are  not  explicitly  compen¬ 
sated  by  the  closed  loop  control  law.  The  extent  to  which  online  implementation  of  PLF 
compensation  can  reduce  nonlinear  coupling  will  be  clearly  evident  by  comparing  the  closed 
loop  responses  of  multiaxis  attitude  control  loops  with  and  without  PLF  compensation.  The 
goal  is  to  provide  distinct  qualitative  changes  in  dynamic  response  representing  nonlinear 
coupling.  Since  mode  locking  is  demonstrated  by  driving  the  system  with  periodic  signals 
and  observing  periodic  responses  one  can  focus  the  experiment  on  tracking  behavior  in  spe¬ 
cific  frequency  bands.  This  permits  considerable  flexibility  in  investigating  the  interaction 
with  internal  nonlinear  resonances  due  to  multibody  and  flexible  mode  interactions.  We  an¬ 
ticipate  that  this  flexibility  will  be  extremely  important  the  interaction  between  the  choice 
of  reduced  models  for  structural  flexure  and  the  design  of  multiloop  tracking  systems. 
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7  Conclusions  and  Directions 

The  principal  conclusions  from  the  three  year  study  are  that  feedback  linearization  offers  a 
comprehensive  approach  to  the  design  and  implementation  of  precision  control  laws  for  rapid 
slewing  and  precision  pointing  of  multibody  systems  with  structural  flexure.  The  methodol¬ 
ogy  for  design  and  implementation  of  control  laws  for  nonlinear  systems  offers  a  fundamental 
approach  which  can  be  readily  employed  for  a  wide  variety  of  system  designs.  Specific  con¬ 
cept  designs  were  completed  and  simulation  studies  revealed  that  robust  control  performance 
can  be  obtained  with  reduced  order  models.  Parametric  uncertainty  in  models  can  affect 
performance  of  PLF  compensation.  In  this  study  adaptive  control  techniques  have  been  de¬ 
scribed  which  can  correct  for  parameter  uncertainty.  We  have  shown  that  certain  standard 
parameter  adaptive  methods  which  have  been  developed  for  linear  systems  (e.g.  model  ref¬ 
erence  schemes)  may  be  inadequate  or  awkward  for  adaptive  PLF  compensation.  For  such 
problems  we  have  shown  an  alternate  adaptive  control  method  which  achieves  asymptotic 
stability  under  much  less  restrictive  assumptions  on  the  structure  of  the  nonlinearities. 

It  is  our  recomendation  that  validation  of  achievable  levels  of  performance  with  PLF  com¬ 
pensation  can  best  be  established  from  laboratory  experiment  employing  actuators  charac¬ 
teristic  of  space  applications.  Providing  for  multiaxis,  large  angle  motions  in  the  experiment 
can  demonstrate  significant  nonlinear  couplings.  We  have  outlined  a  protocol  for  validating 
the  significance  of  these  nonlinear  couplings  and  the  extent  to  which  they  can  be  compen¬ 
sated  using  PLF.  The  protocol  is  involves  observing  steady  state  controlled  motions  with 
periodic  reference  signals  as  inputs  to  the  control  system. 
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